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INTRODUCTION 


In  the  past,  numerous  highway  bridges  have  suffered  extensive 
damages  due  to  strong  motion  earthquakes  [1,2]*.   The  older  bridges, 
consisting  of  single  or  multiple  simple  truss  or  girder  spans  sup- 
ported on  massive  piers  and  abutments,  were  particularly  vulnerable 
to  the  action  of  strong  ground  motions.   Seismic  damages  were  most 
commonly  caused  by  foundation  failures  resulting  from  excessive  ground 
deformation  and/or  loss  of  stability  and  bearing  capacity  of  the 
foundation  soils.  As  a  direct  result,  the  sub-structures  often  tilted, 
settled,  slid,  or  even  overturned;  thus  severe  cracking  or  complete 
failure  was  often  experienced.  These  large  support  displacements  also 
caused  relative  shifting  of  and  damage  to  the  super-structures,  induced 
failures  within  the  bearing  supports,  and  even  caused  spans  to  fall  off 
their  supports.   It  is  significant  to  note  that  very  little  damage 
occurred  to  these  older  structures  as  a  direct  result  of  structural 
vibration  effects. 

Certain  types  of  modern  highway  bridges  may,  on  the  other  hand,  be 
quite  susceptible  to  damage  from  strong  ground  vibration  effects.   This 
fact  became  very  evident  during  the  San  Fernando  earthquake  of  February 
9,  1971,  when  numerous  reinforced  concrete  highway  bridges  suffered 
severe  damages  [3,4].   These  structures  generally  fall  into  two  classifi- 
cations (1)  long,  multiple-span  bridges,  and  (2)  short,  single  or 
multiple  span  bridges. 

The  long,  multiple-span  bridges  usually  consist  of  either  straight 
or  curved  continuous  reinforced  concrete  box  girder  decks  supported  on 
reinforced  concrete  columns  and  diaphragm  abutments.   The  decks  which  are 
continuous  over  many  spans  are  provided  with  intermediate  expansion  joints 


* 
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having  restrainer  bars  tieing  adjacent  spans  together.  Although  the 
overall  bridge  system  in  this  case  is  an  integrated  space  frame, 
significant  structural  discontinuities  result  from  the  presence  of 
these  joints.  The  lateral  rigidity  of  the  overall  continuous  deck 
system  is  therefore  greatly  reduced.  Consequently,  the  lateral  resist- 
ance of  the  total  bridge  structure  is  provided  mainly  by  its  support- 
ing columns.   These  columns  are  often  of  various  lengths  consistent  with 
the  ground  surface  profile.   Structure- foundation  interaction  effects  may 
be  relatively  small  for  the  taller  bridges  but  can  be  quite  large  for  the 
shorter  bridges  in  which  case  they  must  be  included  in  the  analysis.   In 
any  case,  the  dynamic  response  of  abutment  backfills  have  negligible 
effect  on  the  overall  dynamic  response  of  these  bridges. 

During  a  severe  earthquake,  long  multiple-span  bridge  structures  of 
this  type  can  develop  large  amplitude  oscillations.   These  oscillations 
may  cause  large  cyclic  inelastic  deformations,  of  a  coupled  form,  to 
develop  in  the  columns.  Also,  cyclic  slippage  of  the  Coulomb  type  can  take 
place  in  the  expansion  joints  causing  multiple  impacts  and  separations  to 
occur.   These  separations  may  be  sufficiently  large  to  cause  yielding  in 
the  tension  restrainer  bars  and,  as  experienced  in  San  Fernando,  may  even 
cause  deck  spans  to  fall  off  their  supports.  While  yielding  of  the  restrain- 
er bars  can  be  tolerated,  collapse  of  the  bridge  system  is  unacceptable. 

From  the  above  description,  it  is  quite  apparent  that  the  dynamic 
response  characteristics  of  long  multiple-span  reinforced  concrete  bridge 
structures  are  very  complex  due  to  the  nonlinear,  discontinuous  behavior 
in  the  expansion  joints  under  combined  loading  conditions.  Because  these 
characteristics  are  unique,  one  cannot  compare  the  seismic  response  of  this 
type  of  structure  directly  with  the  dynamic  response  of  buildings  under 
similar  excitation  conditions.   Therefore,  it  is  essential  that  the  dynamic 
^response  and  failure  characteristics  of  these  bridges  be  investigated 
separately  by  analytical  and  experimental  means  and  that  the  results  of 
these  investigations  be  fully  correlated  with  field  evidence. 

Turning  our  attention  now  to  the  second  class  of  bridge  structures 
mentioned  above,  i.e.,  short,  single  or  multiple  span  bridges,  it  is 
immediately  apparent  that  structure-foundation  interaction  effects 
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greatly  influence  the  seismic  response  characteristics  of  this  class 
of  structure.   Because  they  lack  intermediate  expansion  joints,  the 
deck  responds  essentially  as  a  rigid  body  transferring  its  seismic  loads 
through  the  supporting  columns  and  abutments  to  the  foundation. 

The  seismic  response  of  abutment  backfills  may  be  in  phase  with 
the  response  of  the  bridge  deck  in  such  a  manner  that  large  dynamic 
active  pressures  develop  which  are  additive  to  the  seismic  forces  in 
the  deck,  thus  greatly  amplifying  the  dynamic  response  of  the  complete 
structural  system.   If  the  bridge  should  be  skewed,  as  is  often  the 
case,  the  active  forces  on  abutment  can  develop  large  twisting  moments 
about  a  vertical  axis  through  the  elastic  center  of  the  bridge.  These 
moments  may  become  sufficiently  large  to  cause  failures  of  columns  and 
abutments  [1,2,3]. 

Obviously,  any  seismic  analysis  of  this  type  of  structure  must 
consider  foundation  and  abutment  backfill  soils  as  part  of  the  complete 
structural  system,  if  realistic  results  are  to  be  obtained. 

As  previously  mentioned,  the  1971  San  Fernando  earthquake  caused 
heavy  damages  to  both  classes  of  modern  reinforced  concrete  bridges 
described  above  which  were  designed  by  the  tradiational  elastic  approach, 
i.e.,  by  the  equivalent  static  seismic  coefficient  method  similar  to 
that  previously  adopted  for  buildings  [5].   Because  of  this  experience, 
it  was  immediately  apparent  that  the  seismic  requirements  used  were 
inadequate  and  that  the  design  methodology  should  be  critically  examined. 
Action  was  quickly  taken  following  the  earthquake  to  correct  certain 
design  deficiencies  [6] ;  however,  the  basic  static  approach  to  design 
still  remains  in  effect. 

Recognizing  the  urgent  need  for  both  theoretical  and  experimental 
research  related  directly  to  seismic  effects  on  bridge  structures,  a  three- 
year  investigation  entitled  "An  Investigation  Of  The  Effectiveness  of 
Existing  Bridge  Design  Methodology  in  Providing  Adequate  Structural  Re- 
sistance to  Seismic  Disturbances"  was  initiated  in  1971  within  the  Earth- 
quake Engineering  Research  Center,  University  of  California,  Berkeley, 
under  the  sponsorship  of  the  U.  S.  Department  of  Transportation,  Federal 
Highway  Administration.   This  investigation  consists  of  the  following 
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five  phases: 

(1)  A  thorough  review  of  the  world's  literature  on  seismic 
effects  on  highway  bridge  structures  including  damages 

to  bridges  during  the  San  Fernando  earthquake  of  February 
9,  1971. 

(2)  An  analytical  investigation  of  the  dynamic  response  of  long, 
multiple-span,  highway  overcrossings  of  the  type  which 
suffered  heavy  damages  during  the  1971  San  Fernando 
earthquake . 

(3)  An  analytical  investigation  of  the  dynamic  response  of 
short,  single  and  multiple  span,  highway  overcrossings 
of  the  type  which  suffered  heavy  damages  during  the  1971 
San  Fernando  earthquake. 

(4)  Detailed  model  experiments  on  a  shaking  table  to  provide 
dynamic  response  data  similar  to  prototype  behavior  which 
can  be  used  to  verify  the  validity  of  theoretical  response 
predictions . 

(5)  A  thorough  comparison  of  dynamic  response  obtained  from 
analyses,  experiments,  and  field  experience  followed  by  the 
preparation  of  recommendations  for  changes  in  seismic  design 
specifications  and  methodology  as  necessary  to  provide 
adequate  protection  against  future  earthquakes . 

The  final  report  covering  Phase  1  has  recently  been  published  [1,2] 
and  the  present  report  is  the  final  report  covering  Phase  2. 

The  primary  objectives  of  the  investigation  in  Phase  2  are  to 
develop  suitable  analytical  procedures  for  linear  and  nonlinear  three- 
dimensional  earthquake  response  analysis  of  long,  multiple- span ,  modern 
highway  bridges  of  the  type  which  suffered  heavy  damages  during  the 
1971  San  Fernando  earthquake  and  to  identify  and  investigate  the  important 
parameters  which  significantly  affect  the  response  of  this  type  of  bridge 
system. 

To  achieve  these  objectives,  linear  and  nonlinear  mathematical 
models  suitable  for  modelling  the  dynamic  characteristics  of  this  type 
of  bridge  have  been  defined  and  analytical  procedures  and  computer 
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programs  have  been  developed  for  determining  the  seismic  response  of 
the  complete  bridge  system  to  arbitrary,  but  prescribed,  earthquake 
excitations. 

Linear  and  nonlinear  seismic  analyses  have  been  performed  on 
three  typical  major  overcrossings  using  earthquake  ground  motions  of 
several  intensities.   These  bridges  are  (1)  the  South  Connector 
Overcrossing  of  the  Golden  State  freeway  and  Antelope  Valley  freeway 
interchange  (5/14  South  Connector  Overcrossing) ,  (2)  the  proposed 
curved  Figueroa  Street  Undercrossing  Connector  in  the  Los  Angeles  area, 
and  (3)  a  straightened  version  of  the  Figueroa  Street  Undercrossing 
Connector.  Various  parameter  studies  have  been  carried  out  for  these 
bridges  and  the  results  are  discussed  and  correlated  with  the  apparent 
behavior  of  this  type  of  structure  during  the  San  Fernando  earthquake . 
In  particular,  the  causes  of  collapse  of  the  5/14  South  Connector 
Overcrossing  are  identified  and  examined. 

Finally,  based  on  the  analytical  results  obtained,  some  general 
conclusions  and  recommendations  related  to  fundamental  seismic  design 
methodology  have  been  formulated. 

Chapter  2  of  this  report  describes  the  linear  and  nonlinear 
mathematical  modelling  of  typical,  long,  multiple-span,  modern  highway 
overcrossings;  Chapter  3  presents  in  detail  a  three-dimensional 
elasto-plastic  flexural  column  model  suitable  for  representing  the 
coupled  yielding  behavior  of  columns  under  large  cyclic  deformations; 
Chapter  4  defines  an  appropriate  nonlinear  mathematical  model  for 
expansion  faults  which  simulate  the  coupled  effects  of  impact  slippage 
with  Coulomb  friction  and  separation  under  the  action  of  elasto- 
plastic  restrainer  bar  forces;  Chapter  5  and  6  describe  the  analytical 
procedures  and  computer  programs;  Chapter  7  presents  the  numerical 
results  of  the  investigation;  Chapter  8  gives  a  discussion  of  the 
numerical  results ,  and  Chapter  9  summarizes  the  conclusions  and 
recommendations . 
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II   MATHEMATICAL  MODELLING  OF  LONG  MULTIPLE-SPAN  BRIDGES 


A.   STRUCTURAL  SYSTEM 

The  structural  system  of  this  type  of  bridge  consists  of  a  multiple- 
span  continuous  box  girder  deck  supported  on  and  rigidly  connected  to 
reinforced  concrete  columns  and  diaphragm  abutments.   The  deck  may  be 
straight  or  curved  in  both  horizontal  and  elevation  views  and  is  sup- 
ported at  discrete  locations  along  its  longitudinal  axis  by  either 
centrally  positioned  columns  or  by  transverse  rows  of  multiple  columns. 
Intermediate  expansion  joints  divide  the  deck  into  several  segments. 
Figure  1  shows  the  structural  lay-out  of  the  5/14  South  Connector 
Overcrossing  which  is  typical  of  this  type  of  bridge. 

The  entire  structural  system  exhibits  characteristics  of  a  continuous 
space  frame.   Its  dynamic  response  to  earthquake  excitations  is  of  lower 
mode  type;  hence,  a  mathematical  model  of  discrete  form  can  be  used  to 
approximate  the  continuous  system.   This  form  of  modelling  leads  to  a 
system  having  a  finite  number  of  degrees  of  freedom. 

Following  the  standard  finite  element  procedure,  these  degrees  of 
freedom  are  chosen  as  the  nodal  displacements  of  the  discrete  finite 
element  model.   For  a  three  dimensional  model,  each  nodal  point  usually 
has  6  degrees  of  freedom,  i.e.,  3  translation  components  and  3  rotation 
components.   Internal  constraints  may  however  reduce  this  number  at  some 
nodal  points. 

For  purposes  of  dynamic  analysis,  the  stiffness,  mass,  and  damping 
properties  of  each  finite  element  must  be  realistically  defined  and  the 
characteristics  of  each  expansion  joint  must  be  properly  modelled.   These 
features  of  the  overall  system  are  described  in  detail  in  the  subsequent 
sections  of  this  report. 

B.   STIFFNESS  IDEALIZATION 

The  finite  element  idealization  of  the  complete  bridge  system 
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results  in  a  stiffness  matrix  which  is  an  assemblage  of  the  generalized 
stiffness  matrices  for  individual  elements,  i.e.,  in  symbolic  form: 

N 

K  -  E  k.  (1) 

—      —i 
1=1 

where  K  is  the  total  stiffness  matrix  for  the  entire  bridge,  k^  is  . 
the  stiffness  matrix  for  element  i,  and  N  is  the  total  number  of  ele- 
ments in  the  system. 

For  small  amplitude  response,  the  bridge  system  may  be  modelled 
by  a  set  of  linear  elastic  elements;  however,  when  subjected  to  high 
amplitude  response  as  occurs  during  severe  earthquakes,  certain 
critical  regions  of  the  structure  may  undergo  large  cyclic  inelastic 
deformations.  Therefore,  nonlinear  finite  elements  must  be  chosen  for 
the  mathematical  model  which  have  realistic  nonlinear  hysteretic  force- 
deformation  characteristics.  The  stiffnesses  of  these  elements  are 
time  dependent  and  are  functions  of  element  deformations  and  deformation 
histories.  Usually,  they  are  linearized  for  analysis  in  a  piecewise 
fashion  using  tangent  stiffnesses  at  discrete  times.  Thus,  the  total 
stiffness  matrix  for  the  entire  structure  may  be  written  in  the  sym- 
bolic form: 

N 
K,.  =  I     K  (2) 

~*   i=i  "fcl 

where  K.  is  the  total  stiffness  matrix  at  time  t  and  k^. .  is  the 
-t  -ti 

stiffness  matrix  for  element  i  at  time  t.   It  should  be  noted  that  all 
nonlinear  elements  referred  to  in  this  report  account  for  material 
nonlinearities  only.   Nonlinearities  arising  from  large  geometry 
changes  are  not  included  as  they  are  negligible  for  the  levels  of  res- 
ponse considered  acceptable  in  these  structural  systems. 

1.  Decks  -  The  decks  of  most  modern  reinforced  concrete  highway 
bridges  are  of  the  curved  box  girder  type  as  shown  in  Fig.  2a.  The 
performance  of  these  girders  under  static  loadings  has  been  investi- 
gated extensively  [7  -  10] .   The  analytical  methods  employed  used 
various  types  of  elements  for  the  deck  structure  such  as  simple 
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beam  elements,  folded  plate  elements,  and  guadralateral  plate  ele- 
ments [7] .   It  has  been  shown  that  when  determining  internal  stress 
distributions  in  constituent  flanges  and  webs  of  box  girders  under 
localized  loadings  that  rather  elaborate  methods  of  analysis,  such  as 
finite  segment  analysis  or  finite  element  analysis,  must  be  used. 
However,  when  the  external  loadings  are  rather  uniformly  distributed 
and  when  only  resultant  forces  on  transverse  cross-sections,  i.e., 
3  components  of  force  and  3  components  of  moment,  are  required,  a 
simple  beam  analysis  is  usually  sufficient  to  yield  accurate  results. 
Thus,  the  mathematical  model  used  herein  for  seismic  analysis  re- 
presents the  deck  by  a  series  of  either  straight  or  curved  beam  elements 
as  shown  in  Fig.  2b. 

The  stiffness  matrix  of  a  linear  elastic  straight  beam  element  is 
found  in  most  modern  textbooks  on  matrix  structural  analysis  [111 
and  the  stiffness  matrix  of  a  linearly  elastic,  circularly  curved, 
beam  element  can  be  found  in  papers  by  Morris  [12]  and  Tezcan,  et.al. 
[13]. 

Since  a  typical  box  girder  deck  is  extremely  stiff  and  strong  in 
comparison  with  its  supporting  columns  and  abutments,  the  high  ampli- 
tude bridge  response  produced  during  severe  ground  shaking  will  be 
caused  primarily  by  deformations  in  the  columns,  abutments,  and 
expansion  joints.  The  deck  will  remain  elastic  and,  therefore,  can 
be  modelled  by  linear  elastic  elements.  Nonlinear  yielding  elements 
must,  however,  be  used  for  columns,  abutments,  and  the  tension 
restrainer  tie  bars  used  in  expansion  joints. 

2.  Columns  and  Abutments  -  The  structural  behavior  of  columns  can 
be  adequately  modelled  using  simple  beam  elements?  see  Fig.  2c. 
Because  of  large  amplitude  response,  coupled  inelastic  flexural  de- 
formations may  occur  in  these  members.  Therefore,  nonlinear  beam  ele- 
ments which  realistically  characterize  the  inelastic  hysteretic  be- 
havior must  be  used.  An  appropriate  elasto-plastic  element  under 
combined  loading  has  been  defined  for  this  purpose  and  described  in 
detail  in  Chapter  3. 

The  abutments  of  this  type  of  bridge  are  usually  of  the  diaphragm 
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type  which  act  in  the  manner  of  shear  wall  elements.   Failures  are 
likely  to  be  of  the  shear  type  causing  excessive  damage.  Therefore, 
such  failures  should  be  avoided  through  proper  design.  Consequently, 
elastic  behavior  is  always  assumed  for  abutments  in  the  present  analysis . 

3.   Foundations  -  Complex  mathematical  models  have  been  used  for 
structure  foundations  in  previous  investigations  [14,  15,  16].   How- 
ever, a  simpler  and  more  approximate  model  has  been  adopted  herein. 
This  model  consisting  of  3  translational  and  3  rotational  springs  is 
used  to  connect  the  base  of  each  column  and  abutment  to  a  rigid  found- 
ation where  the  seismic  excitation  is  fully  prescribed;  see  Fig.  2c. 
For  linear  analysis,  the  stiffnesses  of  these  springs  can  be  evaluated 
using  linear  elastic  half-space  theory  similar  to  that  used  in  the 
reference  [15] ,  or  methods  reported  in  the  literature  [17] . 

For  large  amplitude  response,  the  foundation  soils  may  undergo  in- 
elastic deformations  of  the  hysteretic  type.   In  this  case,  the  six 
foundation  springs  should  be  bilinear  hysteretic  springs.  The  elastic 
stiffnesses  of  these  springs  can  be  obtained  using  standard  methods; 
however,  their  yielding  stiffnesses  and  yield  values  can  be  established 
only  through  extensive  experimental  studies  on  the  dynamic  properties  of 
the  foundation  soils  [15] . 

A  very  simple  method,  which  is  often  used  for  modelling  foundation 
flexibility,  is  to  assume  that  the  columns  and  abutments  extend  below 
the  ground  surface  to  specified  depths  at  which  fully  flexity  conditions 
are  believed  to  have  been  reached.   The  ground  excitations  are  then 
fully  prescribed  at  these  effective  depths. 

C.   MASS  IDEALIZATION 

The  continuous  mass  of  the  bridge  structural  system  is  modelled  in 
discrete  form  by  lumping  element  masses  at  their  end  nodal  points  * 
Since  inertia  forces  are  associated  with  each  of  the  six  degrees  of 
freedom  at  a  nodal  point,  each  lumped  mass  must  be  assigned  appropriate 
moments  of  inertia  about  its  own  three  coordinate  axes  (x,y,z) .   It 
should  also  be  noted  that  when  conducting  nonlinear  dynamic  analyses,  the 
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instantaneous  stiffness  matrix  may  become  singular  in  which  case  it  is 
essential  that  mass  moments  of  inertia  be  assigned  to  each  rotational 
degree  of  freedom.  Following  this  procedure,  a  diagonal  mass  matrix 
m.  is  established  for  each  element  i  (i  =  1,2,  — — ,N) .  The  diagonal 
mass  matrix  for  the  complete  bridge  system  can  then  be  assembled  and 
expressed  in  the  symbolic  forms 

N 
M  =  I     m.  (3) 

i=i  _1 

In  determining  overall  dynamic  response,  this  lumped  mass  method 
has  been  found  to  be  quite  adequate  for  analysis  purposes  [18] . 

D.   DAMPING  IDEALIZATION 

Velocity  dependent  damping  in  the  structural  system  is  represented 
by  a  generalized  damping  matrix  associated  with  the  finite  degrees  of 
freedom  permitted  in  the  mathematical  model.  This  matrix  can  be  derived 
by  consistent  procedures  similar  to  those  used  in  deriving  the  stiffness 
matrix  provided  the  internal  damping  mechanisms  within  each  element  are 
specified.  The  structural  damping  matrix  for  the  complete  structural 
system  would  then  be  of  the  symbolic  form: 

N 
C  =  Z  c.  (4) 

"   i-,  -1 

where  c.  is  the  damping  matrix  for  the  ith  element. 

In  practice,  however,  it  is  very  difficult  to  establish  the  basic 
characteristics  of  damping  in  the  individual  elements.  Therefore ,  it  is 
customary  to  assume  that  the  damping  forces  of  each  element  consists  of 
one  set  which  are  proportional  to  the  velocities  of  each  mass  point  and 
a  set  which  are  proportional  to  the  rates  of  deformation  in  the  element 
[18] .   Thus,  the  element  damping  matrix  can  be  represented  ass 


c.  =  a.  m.  +  $.  k.  (5) 

-i    l  -i    i  -l 
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where  a.  and  8.  are  the  scalar  proportionality  constants  for  the  ith  ele- 
ment.  If  the  entire  structure  is  made  of  the  same  basic  material ,  it 
is  reasonable  to  assume  similar  damping  characteristics  for  each  element 
in  the  system,  i.e.,  a.  and  8.  are  independent  of  i  in  which  case, 
combining  Eqs.  (4)  and  (5),  the  structural  damping  matrix  becomes? 

C  =  a  M  +  3  K  (6) 


where 


a  =  cu  ?  8  -  Si      (i«l,2,— -,N)         (7) 


In  the  case  of  linear  response,  it  can  be  shown  that  the  two 
generalized  parameters  a  and  3  can  be  determined  after  assigning 
damping  ratios  to  the  first  two  natural  modes  of  vibration  [19] . 

For  nonlinear  dynamic  analysis,  the  viscous  damping  properties 
of  the  structure  are  much  more  difficult  to  access.  However,  it  is 
reasonable  to  make  assumptions  similar  to  the  elastic  case,  i.e.,  to 
assert  one  set  of  damping  forces  proportional  to  the  nodal  velocities, 
as  assumed  in  the  elastic  case,  and  a  second  set  which  are  proportional 
to  the  rates  of  elastic  deformation  only.   In  this  case,  the  structural 
damping  matrix  takes  the  form: 


c  =  a  M  +  3  k  (8) 

— t    —    — t 

where  a  and  8  are  the  constants  given  in  Eq.  (7). 


E.   EXPANSION  JOINT  IDEALIZATION 

A  typical  expansion  joint  is  shown  in  Figs.  3a  and  3b.  The  bridge 
deck  on  one  side  of  the  expansion  joint  rests  on  a  ledge  of  the  deck  on 
the  opposite  side.  A  number  of  transverse  shear  keys  exist  on  the  sup- 
port ledge  which  prevent  relative  transverse  motion  of  the  deck  within 
the  joint.  Longitudinal  restrainer  bars  are  attached  across  the  expansion 
joint  tieing  together  the  end  diaphragms  of  the  adjacent  deck  spans. 
These  tie  bars  are  placed  in  such  a  way  that  they  carry  only  tensile 
forces  after  a  small  tie  gap  is  closed  by  joint  separation.  Another 
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small  gap  is  provided  in  the  joint  to  allow  for  thermal  expansion  of 
the  deck  spans.  Vertical  restrainers  are  installed  in  the  expansion 
joint  to  prevent  lifting  of  one  span  off  the  other  along  the  support 
ledge . 

Referring  to  Figs.  3a  and  3b,  the  expansion  joint  characteristics 
may  be  summarized  as  follows? 

(1)  Relative  translation  of  the  deck  spans  along  the  longitudinal 
x  axis  cause  impacts  when  the  gap  between  adjacent  spans  is  closed, 
slippage  on  the  support  ledge  of  the  Coulomb  type,  and  extensions 
of  the  longitudinal  restrainer  bars  when  joint  separations  exceed 
tie  bar  gaps. 

(2)  Relative  translation  of  the  spans  along  the  transverse  y 
axis  is  prevented  by  the  shear  keys. 

(3)  Relative  translation  of  the  spans  along  the  vertical  z  axis 
is  constrained  by  vertical  restrainer  bars. 

(4)  Relative  rotation  of  the  spans  about  the  x  axis  is  constrained 
by  the  vertical  restrainer  bars. 

(5)  Relative  rotation  of  the  spans  about  the  y  axis  is  freely  per- 
mitted. 

(6)  Relative  rotation  of  the  spans  about  the  z  axis  is  coupled  to 
the  relative  translation  of  spans  along  the  x  axis  and  may  be 
partially  or  fully  constrained  by  a  combination  of  closure  of  the 
gap  between  decks,  actions  of  the  longitudinal  restrainer  bars,  and 
the  presence  of  Coulomb  friction. 

Because  of  these  characteristics,  the  expansion  joint  behaves  as  a  non- 
linear discontinuous  element  causing  different  degrees  of  constraint 
during  the  motion  of  the  structure.  Therefore,  the  total  number  of  de- 
grees of  freedom  in  the  complete  system  changes  from  time  to  time  during 
the  response  causing  complexities  to  arise  in  the  nonlinear  dynamic 
analysis  procedures. 

Two  limiting  cases  of  expansion  joint  behavior,  providing  minimum  and 
maximum  constraint,  are  assumed  for  linear  dynamic  analysis.  These  two 
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linear  cases  are  achieved  by  neglecting  the  effects  of  impact  and 
longitudinal  restrainer  bar  constraints  and  by  assuming  either  zero  or 
infinite  Coulomb  friction.  The  infinite  friction  case  prevents  relative 
translation  of  deck  spans  along  the  longitudinal  x  axis  and  prevents 
relative  rotation  about  the  z  axis;  thus,  it  characterizes  the  prototype 
dynamic  behavior  for  low  amplitude  response.   The  zero  friction  case, 
on  the  other  hand,  approximates  prototype  behavior  for  high  amplitude 
response  following  failure  of  the  short  longitudinal  restrainer  bars 
but  prior  to  yielding  in  the  columns. 

The  overall  bridge  structure  is  very  stiff  in  the  infinite  friction 
case  due  to  arching  action  of  the  deck  from  abutment  to  abutment  which 
carries  most  of  the  lateral  loads.   It  is  however,  relatively  flexible 
in  the  zero  friction  case  since  arching  action  cannot  occur  and  the  lateral 
loads  must  be  carried  in  the  supporting  columns. 

F.   LINEAR  AND  NONLINEAR  MATHEMATICAL  MODELS 

In  summary,  the  mathematical  model  used  for  linear  dynamic  analysis 
consists  of  (1)  linear  elastic  straight  beam  elements,  (2)  linear  elastic 
curved  (circular)  beam  elements,  (3)  linear  boundary  spring  elements,  and 

(4)  linear  expansion  joint  elements,  while  the  model  used  for  nonlinear 
analysis  consists  of  (1)  linear  elastic  straight  beam  elements,  (2)  linear 
elastic  curved  beam  elements,  (3)  elasto-plastic  flexural  beam  elements, 

(4)  bilinear  hysteretic  boundary  spring  elements,  and  (5)  nonlinear 
expansion  joint  elements.   Figure  4  shows  the  discrete  parameter  model 
for  the  5/14  South  Connector  Overcrossing. 
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Ill   ELASTO-PLASTIC  FLEXURAL  COLUMN  ELEMENT 
FOR  REINFORCED  CONCRETE  COLUMNS 


A.   BASIC  ASSUMPTIONS 

As  pointed  out  previously,  limited  inelastic  deformations  may 
be  permitted  in  reinforced  concrete  bridge  columns  under  severe  earth- 
quake conditions.   These  deformations  should  however  be  permitted 
only  when  adequate  ductility  and  energy  absorption  capacities  are 
available.   Therefore,  cyclic  yielding  should  be  allowed  in  flexure 
only  where  high  ductility  can  be  provided.   Since  shear  inelastic 
deformations  are  usually  brittle  in  character,  they  should  be  avoided. 
Thus,  it  is  essential  in  providing  high  seismic  resistance  under 
severe  overload  conditions  that  reinforced  concrete  columns  be  sized 
and  detailed  so  that  yielding  will  be  limited  to  flexure. 

Because  reinforced  concrete  is  a  highly  heterogeneous  composite 
material,  the  overall  inelastic  behavior  of  column  elements  should  be 
defined  in  terms  of  generalized  forces  and  their  corresponding  general- 
ized deformations.   Six  generalized  forces  exist  at  any  section  of  an 
element,  namely,  an  axial  force  P,  two  components  of  transverse  shear 

Q  and  Q  ,  torsion  T,  and  bending  moments  M  and  M  ,  Fig.  5.  Axial 
xy     *z  y      z 

normal  stresses  are  caused  by  generalized  forces  P,  M  ,  and  M  while 

y       z 

transverse  shear  stresses  are  caused  by  Q  ,  Q  ,  and  T.   In  constructing 
an  inelastic  element  for  reinforced  concrete  columns,  the  following 
basic  assumptions  have  been  made  (see  Fig.  5)  : 

(1)  The  generalized  force-deformation  relations  of  the  element 
follows  that  of  an  elasto-perfectly  plastic  (elasto-plastic) 
material  model  having  yield  strengths  corresponding  to  the 
ultimate  capacities  of  the  member. 

(2)  The  ultimate  shear  strengths,  Q   and  Q   ,  and  the  ultimate 

xyu      zu 

torsional  strength,  T  ,  are  infinitely  large. 
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(3)  The  ultimate  axial  strength  P  and  the  ultimate  bending 

strengths  M   and  M   are  determined  from  the  axial 
yu      zu 

stress  distribution  present  on  the  cross-section  under 
ultimate  conditions  and  are  independent  of  the  shear 
stresses  caused  by  Q  ,  Q  ,  and  T. 

(4)  The  interaction  among  P  ,  M   ,  and  M   can  be  represented 

u   yu       zu 

by  a  three-dimensional  generalized  yield  surface  in  the 
force  space  spanned  by  P ,  M  ,  and  M  . 

(5)  The  element  retains  linear  elastic  behavior  between  nodal 
points;  however,  elas to-plastic  behavior  as  defined  in 
assumptions  (1)  through  (4)  above  is  permitted  to  occur 
over  the  cross-section  at  each  end  for  a  distance  Ax 
approaching  zero,  i.e.,  the  inelastic  behavior  is  fully 
concentrated  at  the  ends  of  the  element. 

Assumptions  (1)  through  (4)  relate  to  material  properties  within  th 
concentrated  elasto-plastic  hinges  at  the  ends  of  the  elements. 
Between  these  hinges  the  material  remains  elastic.   It  should  be  under- 
stood that  the  above  assumptions  have  been  made  to  approximate  the 
actual  case  in  which  inelastic  deformations  may  occur  at  any  section 
along  a  column.   If  a  sufficient  number  of  elements  are  used,  a  good 
approximation  results.   This  is  due  to  the  fact  that  the  real  case 
is  approached  in  the  limit  with  increasing  numbers  of  elements. 

The  inelastic  column  element  resulting  from  the  above  assumptions 
may  be  constructed  mathematically  by  carrying  out  the  following  three 
steps : 

(1)  Determine  the  ultimate  strengths  P  ,  M   ,  and  M   for  the 

u   yu       zu 

given  section  under  ultimate  conditions . 

(2)  Define  the  generalized  yield  function  which  controls  the 

yielding  interaction  of  P  ,  M   ,  and  M 
1  u   yu       zu 

(3)  Formulate  the  elasto-plastic  tangent  stiffnesses  in  terms 
of  the  element  force  ai*d  corresponding  nodal  displacement 
increments  during  yielding. 
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The  analytical  details  used  in  carrying  out  these  steps ,  as  well 
as'  the  computational  techniques  used  in  the  step-by-step  dynamic 
analysis  procedure,  are  presented  in  the  subsequent  sections  of 
this  chapter. 

B.   ULTIMATE  STRENGTHS 

The  ultimate  strengths  of  a  reinforced  concrete  column  subjected 
to  combined  axial  load  and  biaxial  bending  can  be  determined  using  the 
following  idealizations  [20,21]: 

(1)  Plane  sections  remain  plane. 

(2)  Concrete  carries  no  tension. 

(3)  Maximum  compressive  strain  of  concrete  equals  0.003, 

(4)  The  concrete  stress  distribution  in  the  compression 
zone  can  be  replaced  by  an  equivalent  stress  block 
having  uniform  intensity  equal  to  0.85  times  the  concrete 
cylinder  strength  f '  and  distributed  over  an  equivalent 
compression  zone  bounded  by  the  boundary  of  the  cross- 
section  and  a  straight  line  located  at  a  distance  (1-k  )  c 
from  and  parallel  to  the  neutral  axis  where  c  is  the  dis- 
tance from  the  location  of  maximum  compressive  strain  to 
the  neutral  axis  and  where 


0.85  ,  f '  <  4000  psi 

ki  :    {  /f'-4000\  O) 

0.85-0.05  I -£—— — )  ,  f  >  4000  psi 

\  1000   /  '  c        * 


(5)   The  stress-strain  relation  for  the  steel  reinforcement 
follows  an  elasto-perfectly  plastic  law  having  equal 
yield  stresses  in  tension  and  compression. 
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The  idealized  stress-strain  distribution  at  ultimate  condition 

is  shown  in  Fig.  6  for  a  typical  bridge  column  section.   In  this 

figure,  the  y  axis  is  selected  as  the  minor  axis  of  the  section  and 

the  z  axis  is  selected  as  the  major  axis.   Using  conditions  of 

equilibrium,  the  ultimate  axial  force  P   and  the  ultimate  bending 

moments  M   and  M   can  be  expressed  in  the  form 
yu      zu 


P      =   -0.85    f '    A      -    Z      Af  +      Z      A    f 

u  cc  esse  tsst 


M        =   -0.85    f    A   z      -   Z      Af      z  +      Z      A    f      z   a  ,... 

yu  c      c   c  c      s   sc   sc  t     s   st   st  (10) 

M        =0.85    f '    Ay      +    Z      Af      y  -      Z      A    f      y   j 

zu  c     c   c  c     s   sc  sc  t     s   st  st 


where  the  signs  of  P  ,  M   ,  and  M   agree  with  the  convention  shown 

u   yu       zu 

in  Fig.  6,  the  symbol  Z   sums  over  the  steel  reinforcements  in  the 

compression  zone,  and  Z   sums  over  the  steel  reinforcements  in  the 

tension  zone. 

Since  the  ultimate  strengths  P  ,  M   ,  and  M   of  a  given  column 

u   yu       zu 

section  are  nonlinear  functions  of  the  geometric  parameters  of  the  sect- 
ion, an  iterative  procedure  must  be  used  to  compute  their  values.   An 
initial  position  of  the  neutral  axis  is  first  assumed  in  this  procedure; 
then,  following  an  equilibrium  correction  using  Eqs.  (10),  a  new 
neutral  axis  position  is  determined.   This  procedure  is  repeated  until 
convergence  to  the  correct  neutral  axis  position  is  sufficiently  ob- 
tained. 

Because  modern  bridge  columns  usually  have  rather  complex  shapes 
of  cross-sections  and  the  steel  reinforcements  are  quite  dense  and 
peripherally  arranged,  the  above  described  iterative  procedure  can 
be  very  tedious  if  carried  out  by  traditional  methods.   For  this  reason, 
an  efficient  procedure  using  Newton's  method  was  developed  which  can 
treat  cross-sections  of  arbitrary  shapes  having  arbitrary  distributions 
of  steel  reinforcements.   A  computer  program  YIELD  has  been  developed 
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for  this  purpose 

Having  calculated  the  numerical  values  of  P  ,  M   ,  and  M   for 

U    yU         zu 

a  given  cross-section  and  distribution  of  steel  reinforcements  using 

Eqs .  (10) ,  they  represent  a  point  in  a  three-dimensional  generalized 

force  space  spanned  by  P,  M  ,  and  M  .   If  all  such  possible  points 

were  computed,  they  would  define  a  three-dimensional  surface  which 

controls  the  interactions  among  ultimate  strengths  P  ,  M   ,  and  M 

u   yu       zu. 

It  is  convenient  to  represent  this  surface  by  two  parameters,  namely, 
an  axial  force  ratio  p  and  an  eccentricity  angle  6  as  defined  by 

P 
_  u 

Pu  ~  P 

(11) 

-      6  -  tan"1  (  Ja   ) 

zu 

where  P   is  the  ultimate  axial  compressive  strength  under  conditions 
of  no  bending.   The  interaction  surface  can  then  be  represented  in 
terms  of 

P  =  P   (p  ,6)       ;     P  =  P   (-1,0) 
u    u   u  o    u    ' 

M    =  M    (p  ,0)     ;      M    =  M    (0,0)  (12) 

yu    yu  ^u  yo    yu 

M   =  M    (p  ,6)     ;      M   -  M    (0,5) 
zu    zu   u  zo    zu     2 


For  a  fixed  value  of  p  ,  Eqs.  (12)  define  the  biaxial  bending  inter- 
action curve  between  moments  M   and  M   and,  for  a  .fixed  value  of  9, 

yu      zu 

these  equations  specify  the  interaction  curve  between  the  ultimate 

2       2   I 

axial  force  P  and  the  ultimate  bending  moment  M  =  (M   +  M  )  2 
u  u     yu    zu 

about  an  axis  located  at  an  angle  9  from  the  y  axis.   As  an  example, 
consider  the  typical  bridge  column  section  shown  in  Fig.  7.   Ultimate 
strengths  have  been  calculated  for  this  section  and  interaction  curves 
have  been  established.   Figure  8  presents  a  single  biaxial  bending 
interaction  curve  which  as  shown  is  representative  for  widely  differ- 
ing values  of  p  .   The  ultimate  moments  M   and  M   in  this  case, 
cu  yu      zu 
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as  in  others,  have  been  normalized  with  respect  to  M   and  M   , 

c  yo      zo 

respectively.   Figure  9  shows  a  single  interaction  curve  for  axial 

force  P  and  the  normalized  M   moment  which  is  representative  for 
u  yu  ^ 

6  over  the  full  range  0  <  0  <  90°. 

The  shape  of  the  yield  surface  as  represented  by  Eq.  (10)  is 
obviously  dependent  upon  the  geometric  shape  of  the  cross-section 
and  upon  the  amount  and  distribution  of  reinforcing  steel.   It  is  very 
tedious  however  to  generate  the  complete  yield  surface  for  each  type 
of  section  under  consideration.   Therefore,  for  analytical  purposes, 
approximate  yield  functions  have  been  developed  which  depend  en  fewer 
parameters. 

C.   GENERALIZED  YIELD  FUNCTION 

Many  approximate  yield  surfaces  have  been  proposed  for  reinforced 
concrete  sections  under  combined  loadings  [22,23].   Most  forms  developed 
however  are  for  rectangular  sections  having  simple  arrangements  of  re- 
inforcing bars.   Therefore,  they  cannot  be  readily  applied  to  the  more 
complex  cases  arising  in  modern  bridge  design. 

The  approximate  yield  surface  developed  herein  follows  a  form 
similar  to  that  proposed  by  Bresler  [22]  for  biaxial  bending  and  uses 
a  cubic  polynomial  approximation  for  the  axial  force  and  bending  inter- 
action curves.   As  proposed  by  Bresler,  the  biaxial  bendinq  interaction 
curves  for  fixed  values  of  p  can  be  approximated  by  the  relation 


=  1  (13) 


where  a  and  b  are  constants  which  depend  upon  the  cross-section 

geometries  and  where  M   and  M   are  the  ultimate  bending  moments 

yp      zp 

about  the  y  and  z  axes,  respectively,  for  a  fixed  value  of  p  and 
for  M   and  M   set  equal  to  zero,  respectively.   Moments  M   and  M 

ZU         VU         ^  r  r  i  yP         Zp 

can  be  approximated  by  the  following  cubic  equations: 


M 

a 

M 

_VU 

+ 

zu 

M 

M 

yp 

zp 
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M 

M 
yo 


P  P    2         P    3 

=  1.0  +  ax  (  ^)  +  a2  (  jH)   +  a3  (  ^  ) 
o  o  o 


-  p  <P  <p 
out 


M 


_z£ 


M 


zo 


=  1.0  +\    if)    +b^ (  f)      +b3  (  jS) 
o  o  o 


(14) 


-  P  <P  <P 
out 


where  P   is  the  ultimate  axial  tension  for  the  column  under  conditions 
of  no  bending  and  a,  a  ,  a,  b  ,  b  ,  and  b   are  constants.   Combining 

1     2.  3      *  £  3 

Eqs.  (13)  and  (14) ,  the  approximate  yield  surface  can  be  represented  in 
the  normalized  form 


m 


_£H_ 


1  +  ap  +ap2+ap; 

1*U      2*u       3*u 


m 


zu 


1  +  b.p  +  b0p  2  +  b,p  3 

l^U      2-^u       3^U 


=  1 


or  symbolically 


(15) 


where 


f(p,  ,  m  ,   m  )  =  1 
u   yu   zu 


p=p/P,m   =M/M,m   =M/M 
ru    u'  o'   yu    yu  yo   zu    zu   zo 


(16) 


are  the  normalized  resultant  force  components. 

The  normalized  yield  function,  Eq.  (16),  contains  eleven  constants, 


namely,  P  ,  M 


M 


a,  a. 


,  a3,  b,  bj ,  b 


,   and  b3.   These  constants 


o'   yo'   zo'   '   l'  "2- 

are  determined  by  fitting  smooth  curves  through  ultimate  strength  values 

obtained  by  the  methods  previously  described.   For  typical  bridge  columns 
which  have  nearly  oval-shaped  cross-sections  and  are  peripherally  rein- 
forced, an  ellipse  can  be  used  to  approximate  the  biaxial  bending  inter- 
action curve.   In  this  case,  a  value  of  2  is  used  for  both  a  and  b.   If, 
on  the  other  hand,  the  cross-section  is  a  long,  narrow,  rectangular  section, 
the  corresponding  interaction  curve  approaches  a  straight  line  in  which 
case  a  and  b  both  equal  unity.   Figure  10  shows  a  three-dimensional  plot 
of  the  yield  surface  given  by  Eq.  (15)  which  is  smooth  everywhere  except 
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at  P  =  P  and  P  =  P,_  when  a  and  b  are  both  equal  to  2. 
o         t  ^ 

In  accordance  with  the  definition  of  a  yield  function,  if  the 

values  of  P,  M  ,  and  M  are  such  that  f (p,  m  ,  m  )  <  1,  where 
y      z  y   z 

p  =  P/P  ,  m  =  M  /M   ,  and  m  =  M  /M   ,  then  the  cross-section  is  in  an 
o   y    y  yo       z    z  zo 

elastic  state.   If  p  =  P  ,  m  =  m   ,  and  m  =  m   ,  the  section  is  in  a 

c  u'   y    yu       z    zu 

yield  (or  impending  yield)  state.   When  df  =  o,  the  values  of  p  ,  m   , 
and  m   can  change  but  they  must  always  be  compatible  with  the  yield 
surface.   If  df  <  o,  the  section  is  being  unloaded.   A  set  of  values  for 
p,  m  ,  and  m  which  gives  f (p,  m  ,  m  )  >  1  is,  of  course,  not  possible 
using  an  elasto-perfectly  plastic  model. 


D.   TANGENT  STIFFNESSES  FOR  ELASTO-PLASTIC  FLEXURAL  COLUMN  ELEMENTS 

EP 
In  this  section,  the  elasto-plastic  tangent  stiffness  matrix  k 

which  relates  the  element  force  increments  dS  to  the  element  nodal  dis- 
placement increments  dr  for  an  elasto-plastic  flexural  column  element 
yielding  at  one  or  both  ends  will  be  derived  using  the  generalized  yield 

function  f (p  ,  m   ,  m   )  established  previously  and  the  associated  flow 

ru   yu   zu  *  2 

rule  for  an  elasto-plastic  solid.   The  desired  relation  is 

EP 

dS  =  k*  dr  (17) 

d^l)  [dr-I 

■dsJ  Uj- 

where  dS  ,  dr  ,  K  =  I,  J  are  the  vectors  of  element  force  increments  and 

— K    — K. 

corresponding  element  nodal  displacement  increments  at  node  K. 

By  virtue  of  assumption  (5)  in  Section  A  and  the  assumption  that 
the  element  deformation  increments  can  be  additively  decomposed  into  an 
elastic  component  and  a  plastic  component,  the  element  nodal  displace- 
ment increments  dr  can  be  expressed  as 

dr  =  drE  +  drP  (19) 


where 


-25- 


where 


drE  = 


<£ 


•£ 


drP  = 


< 


< 


(20) 


and  where  dr  ,  K  =  I,  J,  represents  the  vector  of  elastic  nodal  dis- 

"K  P 

placement  increments  at  node  K,  and  dr  ,  K  =  I,  J  represents  the  plastic 

— K 

displacement  increments  resulting  from  concentrated  plastic  deformations 
at  end  K.   The  element  forces  are  uniquely  determined  by  the  element 
elastic  nodal  displacements,  i.e., 


,  E   E 
S  =  k   r 


(21) 


and  therefore 


E     E 
dS  =  k   dr 


(22) 


where  k  is  the  elastic  stiffness  matrix  of  the  element. 

According  to  the  associated  flow  rule  for  an  elasto-plastic  solid 

P 
[24] ,  the  plastic  deformation  increments  dr  at  a  particular  section  K, 

is  governed  by 


<- 


9f_ 

3s 
-KJ 


dX_ 


K 


(23) 


where  dX   is  a  positive  scalar  proportionality  constant  which  can  be 
determined  by  the  yield  condition,  i.e., 


(24) 


f (sK)  -  1 


df 


3f 


-T 


L3§KJ 


dsK=o 


where  T  indicates  the  transpose  operation.   Note  that  if  the  section  is 
not  yielding,  Eq.  (23)  still  applies  provided  9f/3Sj.  is  set  equal  to  0, 
Letting  K  =  I,  J,   the  flow  rule  can  be  expressed  as 
P 


P 
dr  = 


dr*  ^ 

-I 


«z 


if 
3§I 


o 

fdXT  ] 

T 

_ 

"3f  7 

{           } 

ds 

8f 

Kj 

=  *r    dX    (25) 
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and  the  yield  condition  becomes 


f  = 


fCSjM 


f(s-n, 


,   df  = 


9f 
3§I 


3f 


fdS  1 
-I 


id*r. 


9f 
9s 


iT 


dS   =   o 


(26) 


Premultiplying  both  sides  of  the  equation 


E    E     F  P 

dS  =  k  dr  =  k   (dr  -  dr  ) 


(27) 


by  [9f/9s]  ,  one  obtains  upon  using  Eq.  (26)  the  relation 


L  9§ 


-•T 


9f 

9s 


kE   (dr  -  drP)  =  o 


dS  = 
Substituting  Eq.  (25)  into  Eq.  (28),  one  obtains 


(28) 


"  9f  ' 
9S 

T  sE 

r3f  " 

;1  j 

dX     = 

"3f  " 
9s 

kE  dr 


(29) 


Solving  for  dX  gives 


dX  = 


r  9f " 

V 

r  3f  i 

1" 

9f  " 

9S 

kE   dr 


(30) 


which  after  substituting  into  Eq.  (25)  gives 


drP  = 


9f 
9S 


9f 
9s 


9f 
9S 


-l 


9f 
3s" 


kE  dr 


(31) 


Substituting  Eq.  (31)  into  Eq.  (27)  the  required  relation 


dS  =  (kE  -  k*)  dr  =   kEP  dr 


(32) 


-27- 


is  obtained  where 


~t 


raf  ■ 

ft: 

t 
kE 

[&sj 

)-' 

r3f  1 
5J 

(33) 


Physically  speaking,  the  matrix  k  represents  the  losses  of  elastic 

stiffness  in  the  element  due  to  yielding,  and  the  elasto-plastic  tangent 

EP 
stiffness  matrix  k  represents  the  reduced  elastic  stiffnesses  of  the 

EP 
element.   It  should  be  mentioned  that  matrix  k   may  be  singular  depending 

upon  the  location  of  the  point  of  interest  on  the  yield  surface. 


E.   COMPUTATIONAL  TECHNIQUES 

The  evaluation  of  the  elasto-plastic  tangent  stiffness  matrix 
appearing  in  Eq.  (33)  is  straight  forward  theoretically  once  the  yield 
point  representing  the  element  force  state  on  yield  surface  is  deter- 
mined.  However,  in  practice  a  step-by-step  analysis  procedure  is 
usually  required  in  which  case  the  yield  point  determination  is  only 
approximate.   To  be  more  specific,  a  point  inside  the  yield  surface  can 
move  outside  the  surface  during  a  single  time  interval  provided  the 
element  is  in  a  yielding  state. 

To  treat  this  problem  correctly,  it  is  necessary  to  use  an  iteration 
procedure  within  each  time  interval.   This  procedure  however  requires  a 
tremendous  amount  of  computational  effort.   Therefore,  in  the  present 
analysis,  an  approximate  procedure  which  avoids  iteration  is  adopted.   This 
procedure  is  similar  to  one  used  by  Chi  and  Powell  [25] . 

Referring  to  Fig.  11,  assume  that  the  column  element  at  the  beginning 
of  a  time  interval  At  and  at  time  t  is  in  the  elastic  state  as  represented 
by  point  A  inside  the  yield  surface.  Following  the  established  procedures 
a  vector  of  element  nodal  displacement  increments  Ar  is  determined  and  an 
apparent  new  element  force  vector  S   .   is  then  found  using  the  relations 


it+At =  §t +  A^t  -  s-t +  *E  ^t 


(34) 
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and 

f(W  >  1  (35) 

The  method  then  used  to  determine  yield  point  C  proceeds  by  determining 

E 

a  scalar  factor  U, .  o  <  U,  <  1  such  that  ffS,.  +  y,As  )  =  1.   One  can 
l     -   l  — t     l  — t 

then  evaluate  the  elasto-plastic  tangent  stiffness  matrix  for  point 
C.  A  new  element  force  vector  S  .  ,  as  represented  by  point  D,  is 
then  computed  using  the  relation 

=  F  EP 

§t+At  =  S-t   +  ^i  As-t+   (1"V  *t  A^t  (36) 

Usually,  point  D  is  also  located  outside  the  yield  surface  due  to  the 
local  convex  shape  of  the  yield  surface.   A  second  scalar  factor  y 
is  now  introduced  to  bring  point  D  to  point  E  on  the  yield  surface 
along  the  vector  S   .  .   The  final  element  force  vector  is  determined 
by 

§t+At  =  y2  lt+At   '   0<-»2<-1  (37) 

where  y.  is  determined  so  that  f(y  S   »  )  =  1. 

Finally,  at  the  end  of  the  time  interval,  a  new  elasto-plastic 
tangent  stiffness  matrix  is  evaluated,  using  the  new  yielding  force 
vector  S  .  ,  i.e.,  the  vector  to  point  E.  This  new  tangent  stiff- 
ness matrix  is  then  used  as  the  element  stiffness  matrix  during  the 
next  time  interval. 
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Fig.  5  Local  Coordinate  System  For  a  Bridge  Column 
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Fig.  6  Idealized  Ultimate  Stress-Strain  Distribution  For  a 
Typical  Bridge  Column  Section 
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Fig.  7  Cross-Section  of  a  Typical  Bridge  Column 
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Fig.  8  Biaxial  Bending  Interaction  Curve 
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Fig.  9  Axial  Force  and  Bending  Moment  Interaction  Curve 
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Fig.  10  Yield  Surface  of  Typical  Bridge  Column 
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Fig.  11  Determination  of  Yield  Points 
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IV   NONLINEAR  MODEL  FOR  EXPANSION  JOINTS 


A.   IDEALIZATION 

A  nonlinear  mathematical  model  for  simulating  the  dynamic  behavior 
of  expansion  joints  has  been  defined.   This  model  is  idealized  as 
follows  (see  Fig.  12): 

(1)  The  end  diaphragms  of  the  deck  at  each  expansion  joint 
are  rigid. 

(2)  Contact  within  the  expansion  joint  can  develop  at  only 
two  points,  A  and  B,  located  at  the  ends  of  the  rigid 
diaphragms  and  separated  by  a  distance  d  in  the  transverse 
direction  (y  axis) . 

(3)  At  contact  points  A  and  B,  longitudinal  (x  axis)  impact 
springs  having  a  large  stiffness  k  are  attached  to  one  rigid 
diaphragm  leaving  a  small  gap  A  between  the  other  diaphragm 

G 
and  the  springs. 

(4)  Longitudinal  impact  starts  at  points  A  and  B  when  the 
relative  displacements  of  the  two  rigid  diaphragms  close 
the  gap  A_. 

G 

(5)  Relative  slippage  between  the  two  diaphragms  can  take  place 
in  the  longitudinal  direction  at  points  A  and  B  under  the 
action  of  Coulomb  friction. 

(6)  During  periods  of  slippage,  the  Coulomb  friction  forces  at 
points  A  and  B  are  proportional  to  their  respective  vertical 
contact  forces  and  act  in  directions  opposite  to  their  direction! 
of  slippage. 


* 


The  larger  the  stiffness  is,  the  closer  that  true  impaet  eondition  is 
attended.  The  practical  limit  of  this  value  depends  upon  the  number  of 
significant  digit  carried  by  the  computer  used  and  the  time  interval 
chosen  in  the  analysis. 

-37- 


(7)  A  total  of  N  longitudinal  restrainer  ties  having  equal 
elastic  tensile  stiffnesses  k  and  zero  compressive 
stiffnesses  are  placed  through  the  two  rigid  diaphragms 
with  a  tie  gap  A  at  one  end  of  each  bar  and  at  transverse 

distance  y.  (i=l,2, -,N  )  from  the  longitudinal  center 

line  axis  of  the  deck. 

(8)  Longitudinal  restrainer  ties  are  elasto-plastic  in  tension 
and  have  equal  yield  forces,  S  . 

(9)  Relative  transverse  motions  of  the  two  rigid  diaphragms  are 
restrained  by  a  very  stiff  elastic  shear  spring  of  stiffness 
k   located  at  the  center  point  C. 

(10)  At  contact  points  A  and  B,  the  rigid  diaphragms  are  inter- 
connected by  vertical  (z  axis)  springs  having  stiffnesses  equal 

tO         l^. 

(11)  To  provide  generality  in  the  model,  the  expansion  joint  may 
be  skewed  with  respect  to  the  deck,  i.e.,  the  transverse  axis 
of  the  point  (s  axis;  Fig.  13)  may  not  be  normal  to  the  x  axis. 


The  behavior  of  the  idealized  expansion  joint  shown  in  Fig.  12  can 
be  characterized  conveniently  using  displacement  coordinates,  r,  as  shown 
in  Fig.  13 ,  i.e. 


r  = 


r-I 


T-j) 


r-K  = 


"rAx_ 

r 

y 

rAz 

rBx 

e 

s 

_rBz_ 

(K  =  I, J) 


(38) 


where  9   is  a  joint  rotation  about  the  s  axis  which  may  be  skewed  by  an 
angle  i|i  as  shown.   The  expansion  joint  coordinates  r  can  be  related  to 
the  local  nodal  coordinates  r  as  defined  by 
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M 
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Lfl.. 


(K  =  I, J) 


(39) 


through  the  following  transformation: 


r  = 


?ii 


lij 


a    0 


0    a 


r-I 


r-J 


=   A  r 


(40) 


where  transformation  matrix  a  is  given  by 
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(41) 


The  forces  associated  with  the  expansion  joint  coordinates  are 


S  =• 


V 


■  Sj; 


§k  = 


"SAx" 

s 
y 

SAZ 

SBX 

M 

s 

.SBZ. 

,  (K  =  I, J) 


(42) 


while  the  forces  associated  with  the  local  nodal  coordinates  are: 
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s  =. 


S-I 


S-J 


h- 
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M- 


K 


(K  =  I, J) 


(43) 


It  is  easy  to  show  that  S  are  related  to  S  by  the  transformation 


<h 


S  = 


*J 


T 
a     0 


Si 


(.SjJ 


T  - 
A   S 


(44) 


B.   PIECEWISE  LINEAR  EXPANSION  JOINT  STIFFNESSES 

The  idealized  expansion  joint  can  be  characterized  by  a  stiffness 
matrix  which  is  piecewise  linear  with  a  number  of  discontinuities  de- 
pending upon  the  relative  displacements  of  the   joint.   For  a  step-by- 
step  dynamic  analysis,  if  it  is  assumed  that  these  finite  number  of 
discontinuities  occur  at  the  end  of  each  time  step,  then  the  vector  of 

expansion  joint  force  increments  As  can  be  related  to  the  vector  of  the 

—  — EJ 

expansion  joint  displacement  increments  Ar  by  a  stiffness  matrix  k 

which  is  constant  within  each  time  step,  i.e., 


AS  =  k^J  Ar 


(45) 


— EJ  — 

where  k   is  a  function  of  r  at  time  t. 

-EJ  . 
It  is  convenient  to  define  the  stiffness  coefficients  of  k   in 

-t 

terms  of  the  relative  expansion  joint  displacements  u  defined  by 
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(46) 


and  the  corresponding  expansion  joint  forces  F  defined  by 


F  = 


"FAx_ 

_SAx" 
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(47) 


Thus,  the  idealized  expansion  joint  can  be  characterized  by  a  stiffness 
matrix  k  relating  AF  to  Au,  i.e., 


AF  =  ic  Au 


(48) 


Stiffness  matrix  k  represents  contributions  from  two  sources.   One 
contribution  is  the  stiffness  matrix  k  of  the  idealization  expansion 
joint  without  longitudinal  restrainer  ties  and  the  other  contribution 
is  the  matrix  k  due  to  the  presence  of  longitudinal  ties,  i.e., 


k  =  k,  +  k 

-     -1     -2 


(49) 


The  stiffness  coefficients  of  k  can  be  written  in  the  form 
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(50) 


where 


=  k. 


=  0 


for 

for 

or 


u   <  -  A„  and  u,   <  0, 

Ax  G        Ax  - 

u   <  -  A„  and  u,   >  0, 

Ax  G        Ax 

Ax  G 


=  0 


for 

for 

or 


u   <  -  A  and  u   <  0, 

Bx      G  Bx  _ 

u„  <  -  A„  and  u„  >  0, 

Bx      G  Bx 

Bx  -    G 


(51) 


and  where  k  ,  Jc  are  those  quantities  defined  previously. 

The  stiffness  coefficient  of  each  longitudinal  restrainer  tie  within 
a  time  step  can  be  expressed  as 


AFm.  =  km.Aum.    /   (i=l, ,NJ 

Ti    Ti   Ti  T 


(52) 


where  AF   and  Au   are,  respectively,  the  increment  of  tie  force 

F„.  and  the  increment  of  tie  deformation  um.  in  the  ith  tie  during 
Ti  Ti 

each  time  interval.   The  stiffness  k  .  is  the  instantaneous  stiffness 

Ti 

of  the  ith  tie  bar  at  time  t,  which  can  be  expressed  as 
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=  0    for 


u.  <  (A+uf,.) 

Ti     T  Ti     t 


km.   =  km   for    (Am+u*\)  <  um.  <  (A  -Hi^.+uf),     (53) 
Tl        T  T   Tl      Tl      T   Tl   T 
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P    E 
um.  >  (A_+um.+u  ) 

Tl      T   Tl   T 


P  E 

where  u   >  0  is  the  current  plastic  elongation  of  the  ith  tie,  and  u 

is  the  elastic  elongation  of  the  ties,  i.e.,  u  =  S  /k  . 

Since  the  end  diaphragms  are  rigid,  the  tie  deformation  u  .  can 

be  determined  from  the  longitudinal  expansion  joint  deformations  u   and 

u   by  the  relation 

BX 
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,  (i=l,2,  — ,NT) 


Bx 


(54) 


From  conditions  of  equilibrium,  the  longitudinal  expansion  joint  forces 

F,   and  F_  <-an  be  expressed  as  the  sum  of  all  tie  bar  forces  Fm>,  i.e., 
Ax      Bx  r  Ti 
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(55) 


Combining  Eq.  (52),  (53),  and  (54),  the  following  relation  is  obtained 
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Thus,  the  stiffness  matrix  k.  in  Eq.  (44)  can  be  expressed  as 
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— EiJ 
Using  Eqs.  (46)  and  (47),  the  stiffness  matrix  k   can  now  be 

expressed  as 


r 


k  -k 
-k   k 


(59) 


-EJ 
Finally,  the  stiffness  matrix  k   must  be  transformed  to  the 

local  nodal  coordinate  system  using  Eqs.  (40)  and  (44).   This  trans- 
formation results  in  the  relation 


AS  =  k^J  Ar 


.  EJ      T   -EJ 
k   =  A   k   A 


(60) 
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EJ 
Matrix  k   is  the  stiffness  matrix  which  relates  the  vector  of  local 

nodal  force  increments  As  to  the  vector  of  local  nodal  displacement 

increments  Ar  at  time  t,  and  which,  after  local  to  global  coordinate 

transformation,  can  be  used  to  assemble  the  total  stiffness  K  for  the 

complete  structural  system. 

At  the  end  of  each  time  interval  during  the  analysis  process, 

joint  conditions  as  expressed  by  Eq.  (51)  and  the  condition  of  each 

longitudinal  restrainer  tie  are  checked  and  expansion  joint  stiff- 

p 
nesses  are  formed  accordingly.   The  quantities,  u   ,  u   ,  and  u  . , 

(i=l,2, ,N  )  are  carefully  monitored  since  the  maximum  value  of  u 

and  u   represent  the  maximum  expansion  joint  separations  at  points 

P 
A  and  B  and  the  maximum  values  of  u   ,  (i=l,2, ,N  )  indicate  the  maximum 

ductilities  that  the  longitudinal  restrainer  ties  should  provide. 


COULOMB  FRICTION  FORCES 


Coulomb  friction  forces  C   and  C   ,  which  develop  at  contact 
points  A  and  B  when  the  expansion  joint  undergoes  longitudinal 
relative  displacements  and  when  the  vertical  contact  forces  F   and 
F   are  compressive,  can  be  considered  as  pairs  of  self -equilibrating 
forces  acting  on  the  expansion  joint  along  the  longitudinal  x  axis 
at  points  A  and  B.   According  to  the  assumptions  stated  previously, 
these  forces  can  be  expressed  as 


CAx  =  V<FAz>l  FAzl  Si^n  (*AX> 
CBx  =  V<FBz>l  FBzl  Si*n  (*BX} 


Ax     N  Az  /  I   Az  I    y     Ax 

(61) 


where  v  is  a  constant  coefficient  of  friction  and  where 


=   1  f        <   0 

<FAz>  (62) 

^AZ/=0  F>0 

Az   ~ 

=    1  F        <    0 

<    FBZ>       „  „  '«» 

=   0  F„      >    0 

Bz 
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If  it  is  assumed  that  the  relative  velocities  u   and  u  ,  and  the 

forces  F   and  F   do  not  change  signs  during  a  time  interval  At, 
Az      Bz 

then  the  changes  of  Coulomb  friction  forces  during  the  interval  can 
be  expressed  by 


AC 
Ac 


Ax  "   v  Az^>  ' 


Bx 


AF 


Bz 


sign  (Uax) 


sign  (uBx) 


(64) 


Let  AP  represent  the  incremental  Coulomb  friction  force  vector 

in  the  expansion  joint  coordinate  system  during  time  interval  At.   Then 

_C 
AP  can  be  expressed  as 


AP°  = 


Ap°  n 


APC 
-J 


AP  j  = 


AC 
0 
0 
Ac. 

0 
0 


Ax 


'Bx 


€ 


-AC. 

0 

0 
-AC 

0 

0 


Ax 


Bx 


(65) 


Thus,  the  incremental  Coulomb  friction  force  vector  AP  in  the  local 
nodal  coordinate  system  can  be  obtained  using  Eq.  (44);  i.e. 


APC  =  AT  APC 


(66) 


When  AP  for  all  expansion  joints  in  the  bridge  system  are  transformed 
to  the  global  coordinate  system  and  are  assembled,  an  incremented  Coulomb 
friction  force  vector  AR  in  the  global  coordinate  system  is  obtained. 

Q 

The  vectors  AR  are  then  added  to  the  incremental  external  force  vectors 
during  each  time  interval. 
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Fig.  13  Coordinate  Systems  For  Expansion  Joint 
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V   ANALYTICAL  PROCEDURES  FOR  EVALUATION  OF  SEISMIC  RESPONSES 

A.   EQUATIONS  OF  MOTION 

The  equations  of  motion  for  an  n  degree  of  freedom  system 
expressing  dynamic  equilibrium  at  time  t  can  be  expressed  in  the 
standard  matrix  form 


M  Ufc  -+-  CtU   +  Ktufc  =  R(t)  (67) 

where  M,  C  ,  and  K  are  the  mass,  damping,  and  stiffness  matrices, 
respectively,  and  where  R(t)  is  the  applied  dynamic  load  vectors. 

Vectors  u.  ,  u.  ,  and  uv  are  the  absolute  acceleration,  velocity,  and 

-t  ->t      -t  ■* 

displacement  vectors,  respectively,  as  measured  with  respect  to  a  fixed 
set  of  global  coordinate  axes. 

If  the  system  is  subjected  to  prescribed  support  excitations, 

Q 

a  complete  set  of  nodal  displacements  u  should  be  considered  which 
include,  in  addition  to  the  n  free  nodal  displacements,  the  n  prescribed 
non-zero  support  displacements.  Thus,  the  complete  nodal  displacement 
vector  can  be  expressed  in  the  partitioned  form 


(68) 


where  u  is  a  vector  containing  the  n  non-zero  support  displacements 
(translations  and/or  rotations) .  Vector  u  can  be  conveniently  de- 
composed  into  a  quasi-static  displacement  vector  u  and  a  dynamic  dis- 
placement vector  u  ,  i.e. 


(69) 


where  by  definition  u  =0.   Enlarging  the  mass,  damping,  and  stiffness 

matrices  as  well  as  the  dynamic  load  vector  in  Eq.  (67)  to  account  for 

b 
the  n  support  displacements,  the  equations  of  motion  for  the  complete 
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system  can  be  expressed  in  the  partitioned  matrix  form 


M 


b  T 

(mV 


M* 


M 


bb 


2t 


..b 
St 


St 


<c5T 
-t 


St 

<"b-> 


-t   i   -t 


«*>T!  k»> 

t   i   -t 


u 
-t 


D  / 
U  I 
~t 


~RB<t> 


(70) 


The  equations  of  motion  associated  with  the  n  free  nodal  displacements 
now  become 


'V 


M  --,♦  [^kV  MJk--  » 


St 


R(t) 


u 
-t 


.b 


b 
Ht 


(71) 


Substituting  Eq.  (69)  into  Eq.  (71),  one  obtains 


u 

-s 


5?  H  +  2t  H  +  5 1  H  =  5(t) 


-  [■  ^]  J"  k  £] 


-b 

u 

-s 


u 
-s 


u 
-s 


(72) 


after  making  use  of  the  relation 


b  b 
K.  u  +  K^  u  =0 
-t  -s   -t  -s   -  r 


(73) 


which  is  satisfied  by  definition  of  the  quasi-static  vector  u  .   This 
vector  can  thus  be  obtained  directly  from  Eq.  (73),  i.e., 


-i  Jo     b  _  _ 
u  =  -K.   K~  u 


-s 


-B.  u 
t   -t  -s    -t  -s 


(74) 


where  B  =  K   If  is  a  matrix  of  quasi-static  influence  coefficients 
resulting  from  the  n  non-zero  support  displacements.   If  the  system 


-50- 


is  linear,  all  coefficients  in  matrix  B  are  invariant  with  time. 
Usually  the  damping  terms  on  the  right  hand  side  of  Eq.  (72) 
are  small  compared  to  the  inertia  terms  and  therefore  may  be  dropped 
from  the  equation  without  introducing  significant  errors.  Also,  the 
coefficients  in  matrix  M  can  be  set  equal  to  zero  since  mass  coupling 
vanishes  for  a  lumped  mass  model.   Thus,  after  substituting  Eq.  (74) 
into  Eq.  (72),  the  equations  of  motion  reduce  to  the  form 

H  li  +  C,  u  +  K  u  =  R(t)  +MB.U  (75) 

-  -   -t  -   — t  -   -      —  — t  -s 

where  u  is  a  vector  containing  the  prescribed  support  excitations. 

1.  Multiple  Support  Excitations  -  When  ground  excitations 

corresponding  to  each  of  the  n  support  displacements  are  prescribed 

by  a  vector  u  (t) ,  the  vector  u  in  Eq.  (75)  can  be  expressed  as 
g  — s 

u  =  u  (t)  (76) 

-s   -g 

and  the  equations  of  motion  become 

M  u  +  C.  u  +  K.  u  =  R(t)  +MB,  um(t)  (77) 

-  -    -t  -    -t  -    -       -  -t  -g 

2.  Rigid  Support  Excitations  -  When  ground  excitations  at  all 

support  points  along  the  base  of  the  structure  are  identical  and  are 

••r 

prescribed  by  a  rigid  ground  acceleration  vector  u  consisting  of 

three  translational  components  u  ,  u  ,  and  u  ,  measured  along 
their  corresponding  global  axes  X,  Y,  and  Z,  i.e., 

V 
fg  -  \  > , 

the  equations  of  motion  become 


M  H  +  2t  H  +  £.  u  =  R(t)  +  M  Br  ur  (t)  (79) 
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reliable  phase  difference  predictions  can  be  made.   Therefore,  ground 
motions  in  the  form  of  rigid  base  motions  are  used  in  analyses. 

C.   LINEAR  ANALYSIS  PROCEDURES  -  MODE  SUPERPOSITION  METHOD 


When  the  structural  system  is  linear,  the  stiffness  and  damping 
ces  are  invariant  with  time, 
Eq.  (81)  takes  on  the  linear  form 


matrices  are  invariant  with  time,  i.e.,  K  =  K  and  C  =  C.   Thus 


Mu  +  Cu  +  Ku  =  R(t)  +  M  B  u  (t)  (82) 

A  standard  procedure  can  be  followed  in  obtaining  dynamic  response, 
namely,  solving  the  generalized  eigenvalue  problem  for  mode  shapes 
and  frequencies,  solving  a  decoupled  set  of  normal  equations  of 
motion,  and  using  mode  superposition  to  obtain  time  histories  of 
response  [18] . 

1.   Mode  Shapes  and  Frequencies  -  The  desired  undamped  free  vibration 
mode  shapes  and  corresponding  frequencies  can  be  obtained  by  solving 
the  equation 

K  $.  =  OjJ  M  $.  (i=I,2, ,q);  q  <  n  (83) 

where  CO.  and  $.  are  the  frequency  and  shape  vector,  respectively, 
for  the  ith  mode  and  where  q  is  the  number  of  lowest  modes 
required  for  the  accuracy  of  solution  desired.   Usually,  q  is 
much  less  than  n  and  can  be  efficiently  determined  using  a 
determinant  search  technique  combined  with  inverse  iteration 
or  a  subspace  iteration  technique  as  recently  developed  by 
Bathe  [26] . 

The  modal  matrix  $  =  [^  $,---  $  1  must  satisfy  the  ortho- 
gonality condition 

T        2 
*  ■  K  *  -  0  (84) 
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where  fi  is  a  diagonal  matrix  containing  the  squared  frequencies 

W2,  to2, ,  co2.   It  is  convenient  to  normalize  the  modal  matrix 

.12       q 

so  that  it  satisfies  the  condition 


T 
$  M  $  =  I  (85) 


where  I  is  the  unit  (or  identity)  matrix. 

Vectors  u(t),  u(t),  and  u(t)  can  be  expressed  in  terms  of  the 
modal  matrix  and  the  normal  coordinate  (mode  amplitude)  vector  U(t) 
as  follows: 

u=$\J;  u  =  $U;  u  =  $\3  (86) 

Substituting  Eq.  (86)  into  Eq.  (82) ,  premultiplying  the 

T 
resulting  equation  by  $  ,  and  making  use  of  Eq. (84)  and  (85) , 

one  obtains 

U  +  A  u  +  Q2   U  =  R*(t)  (87) 

where  A  is  a  diagonal  matrix  containing  the  terms  2£.to., 

for   i=l,2. ,q  and  where  R* (t)  is  a  normal  load  vector 

defined  by 


(88) 


R*(t)  =  $T  [  R(t)  +  M  B  u  (t) 
Using  Eq.  (6),  matrix  A  can  be  expressed  in  the  form 

A  =  $T  c  $  =  $T[  CtM  +  $K  ]$  =  ai  +  Bft2  (89) 


2.   Response  Time  Histories  -  The  solution  of  Eq.  (87)  can 
be  carried  out . in  the  time  domain  using  the  convolution  integral 
[18] 

.    t   ,      -5  W  (t-T) 

U.  (t)  =~      /  R.  (T)  e  sin  ol.  (t-T)  dx      ^90) 

l       li)_  .  n   l  Dl 

Di  ° 

for  i=l,2, ,q; 
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where  CO  .  is  the  damped  frequency 

l 

0>  .  =  u.  (1-C2)2  (91) 

Dl     l     l 

Thus,  the  response  time  histories  u(t),  u(t),  and  u(t)  can  now 
be  obtained  using  Eq.  (86) . 


3 .   Spectral  Responses  -  For  a  given  single  component  of  earth- 
quake ground  acceleration  u  (t) ,  the  displacement  response  spectra 

S,  (oj,£)  define  the  maximum  (or  extreme)  values  of  displacement 
a 

governed  by  the  equation  [27] 

V  +  2^   0  +  CO2  U  =  u  (t)  (92) 

g 

These  spectra  are  usually  obtained  directly  from  the  pseudo-velocity 
response  spectra  as  defined  by 

S  (0),?) 
S,  (G0,£)  =   -^— (93) 

u  0) 

Suppose  these  response  spectra  are  known  for  each  ground  acceleration 
input  in  the  vector  ii  (t)  appearing  in  Eq.  (88)  and  suppose  the 
R(t)  =  0/  then  the  modal  maximum  (or  spectral)  response  can  be  ob- 
tained from  the  relation 

|  U.  (t)  |      =  <£   M  B  S,  (0J.,£.)     (i=l,2, ,q)     (94) 

1   i     'max.    ^l  -  -  -d  i   l  ^ 

The  maximum  displacement  response  |u. (t) |     can  then  be  determined 

X      IT13.X  • 

using 

I  u.  (t)  |      =  <J>.  |u.(t)|  (i=l,2,  —  ,q)     (95) 

-l     'max.   -i  '  l    'max.  ^ 

Thus,  the  root-mean- square  spectral  displacement  response  can  be 
expressed  as 
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u(t) 


R.M.S, 


q      .2 
1      ui  It) 

.._...  '  -i   I  max, 


(96) 


D.   NONLINEAR  ANALYSIS  PROCEDURE  -  STEP-BY-STEP  INTEGRATION 

When  the  structural  system  is  nonlinear,  the  coupled  equations 
of  motion,  Eq.  (81) ,  must  be  solved  by  step-by-step  integration 
methods  since  mode  superposition  no  longer  applies.   To  carry  out 
this  integration,  the  equations  of  motion  have  been  formulated 
on  an  incremental  basis.   Considering  a  time  interval  At  starting 
at  time  t  and  assuming  the  stiffness  and  damping  matrices  at 
time  t,  i.e.,  K  and  C  ,  can  be  applied  over  the  full  time 
interval ,  one  obtains  the  equations  of  motion  in  the  incremental 
form 


M  Au  +  C,  Aii  +  K  Au  =  A  R(t)  +  M  B  Au  (t) 
_   _    _t   -    -t   -      -      -g 


(97) 


where 


Au  =  Au(t)  =  u(t+At)  -  u(t) 
Au  =  Au(t)  =  u(t+At)  -  u(t) 
Au  =  Au(t)  =  u(t+At)  -  u(t) 


(98) 


and  where 


AR(t)   =  R(t+At)  -  R(t) 

Aii  (t)  =  u  (t+At)  -  u  (t) 
-9     -g        -g  ' 


(99) 


To  carry  out  the  step-by-step  numerical  integration  of  Eq.  (97) ,  it 
is  necessary  to  introduce  an  approximate  operator  in  the  time  domain 
to  replace  the  differential  operator.   In  the  present  analysis, 
the  Newmark  generalized  acceleration  method  [28]  is  adopted  which 
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assumes  the  following  approximations  for  nodal  velocities  and 
displacements : 

u(t+At)  =  u(t)  +  [(1-6)  u(t)  +  6  u(t+At)] 

(100) 
u(t+At)  =  u(t)  +  u(t)  At  +  [  (|  -  a)u(t)  +  a  u(t+At)] 

where  parameters  6  and  O   can  be  chosen  to  give  the  required  inte- 
gration  stability  and  accuracy.   When  0=2  an(^  a  =  ?'  t*ie 
approximations  correspond  to  the  linear  acceleration  method,  and  when 
6  =  -  and  O   =  r  they  correspond  to  the  constant  acceleration 
method.   While  the  linear  acceleration  method  is  conditionally 
stable  depending  upon  the  magnitude  of  At,  the  constant  acceleration 
method  is  unconditionally  stable  for  any  magnitude  of  At  [29,30]. 
Using  Eq.  (97) ,  the  approximations  given  by  Eq.  (100)  can  be 
expressed  in  the  incremental  form 


Au(t)  =  Cj  Au(t)  -  C3  u(t)  -  C^  u(t) 
Au(t)  =  C2  Au(t)  -  C^  u(t)  -  C5  u(t) 


(101) 


where 


c  =  A      '     c  =  -2-     '     c  =  4- 

1   At2  c2   At  3   At 


C,  =  2  ,  C5  =  0 


(102) 


for  the  constant  acceleration  method  and  where 


c     _  § ,       t  r      _   3_       ,       r  _  6_ 

1   At2  u2  "  At  a  "  At 

C  =  3  C  =  £■ 


(103) 


for  the  linear  acceleration  method. 

Substituting  Eq.  (101)  into  Eq.  (97) ,  one  obtains  after  some 
manipulation 
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C,  M  +  C2  Ct  +  Kt 


Au  =  AR(t) 


(104) 


where 


AR(t)  =  AR(t)  +  M  B  AU  (t)  +  icM  +  CC  )  u(t) 
-g      13-    tf-tj  - 

+  (C.M  +  C5Ct}u(t) 


(105) 


Using  Eq.  (8),  i.e.,  C   =  aM  +  6K  ,  and  introducing  the  following 
normalized  constants 


C?  =  1/(1  +  6  c2) 


c9  =  c7  (c3  +  ac,) 


Cll=C7  (C,  +aC5} 


C8  =  C?  (C,  +  aC?) 


C10=  ^7 


ci2=  ec5c7 


(106) 


Eq.  (104)  can  be  put  into  the  form 


|  C8M  +  KtJ  Au  =  C7  j^AR(t)  +  M  B  Au  (t)J 


(107) 


Defining  an  effective  displacement  vector  Au(t)  by  the  relation 


Au(t)  =  Au(t)  -  C1Q  u(t)  -  C12  u(t) 


(108) 


and  defining  two  new  constants 


C13  E  C9  -C8  C10 


C    3C..-C  C ,  .„ 

1  «♦      11     8    12 


(109) 


Eq. (107)  can  be  written  in  the  form 


K  Au  =  AR(t) 


(110) 
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where 


K   =  C„  M  +  K,  (111) 

-t     8  -    -t 


and 


AR(t)  =  C,  \   AR(t)  +  M  B  u  (t)|  +  C,  ,  M  u(t)  +  C,  ,  M  u(t) 
—         7  (_   -        —  —  — g    J       13—  —  14-  — 

(112) 

Matrix  K  is  the  effective  dynamic  stiffness  matrix  and  AR(t) 
is  the  effective  load  vector. 

Equation  (110)  can  be  solved  for  Au(t)  at  each  time  instant 
t  and  Eqs.  (108)  and  (101)  can  be  used  to  obtain  Au(t),  Au(t), 
and  Aii(t).   The  displacements,  velocities,  and  accelerations  at 
time  t  +  At  can  then  be  determined  from  Eq.  (98) . 

The  displacement  solutions  u(t  +  At)  can  be  used  to  calculate 

the  internal  force  vector  S(t  +  At)  and  the  new  tangent  stiffness 

matrix  k     .   for  each  nonlinear  element  in  the   system.   The  new 
-t  +  At  J 

total  tangent  stiffness  matrix  K    .   is  then  obtained  by  the 
standard  assemblage  procedure.   It  is  important  to  note  that 
when  determining  the  element  internal  forces,  the  static  forces 
existing  in  the  element  must  always  be  included  since  the 
calculation  of  element  tangent  stiffnesses  depend  upon  the 
magnitudes  of  the  element  forces. 

The  step-by-step  integration  algorithm  for  the  nonlinear 
system  presented  above  may  be  summarized  as  follows: 

(1)   Initial  calculations 

(a)  Form  the  initial  stiffness  matrix  K  and  the  mass 
matrix  M  for  the  system. 

(b)  Solve  for  the  initial  displacements  and  the  element 
forces  due  to  static  loads. 

(c)  Compute  the  initial  quasi-static  influence  matrix 
B,  if  required. 
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(d)  Set  up  the  dynamic  load  and  ground  excitation 
time  histories  R(t)  and  u  (t) . 

-g 

(e)  Calculate  the  step-by-step  integration  constants, 
C^  i=l,2, ,14. 

(2)   For  each  time  increment  At,  at  time  t 

(a)  Form  the  tangent  stiffness  matrix  K  m 

(b)  Compute  the  current  quasi-static  influence  matrix  B  , 
if  required. 

(c)  Form  the  effective  dynamic  stiffness  matrix  K  , 

(d)  Triangularize  K  ., 

(e)  Form  the  effective  load  vector  AR(t), 

(f)  Solve  for  the  effective  displacement  increments  Au(t) 

(g)  Compute  current  accelerations,  velocities,  and  dis- 
placements u(t),  u(t),  u(t). 

(h)   Calculate  current  element  forces,  check  nonlinearity 
conditions,  and  compute  new  element  tangent  stiffness 
matrices,  inelastic  deformation  vectors,  and  unbalanced 
force  vectors,  if  necessary. 

(i)   Return  to  step  (a)  or  (e) ,  as  necessary,  for  next 
time  increment. 


-60- 


VI   LINEAR  AND  NONLINEAR  SEISMIC  ANALYSIS  COMPUTER  PROGRAMS 

A.   LINEAR  ANALYSIS  COMPUTER  PROGRAM  -  BSAP 

A  computer  program  BSAP  has  been  developed  for  dynamic 
response  analysis  of  linear  bridge  structural  systems.   This  pro- 
gram is  a  modified  version  of  the  SAP  program  originally  developed 
by  Wilson  and  co-workers  for  static  and  dynamic  analyses  of  linear 
structural  systems  [31] .   Although  the  original  version  of  SAP 
is  very  efficient  and  can  be  used  for  analyzing  a  wide  variety 
of  structures,  some  modifications  are  necessary  to  make  it 
suitable  for  the  analysis  of  bridge  structural  systems.   While 
most  of  BSAP  follows  the  organization  of  SAP  [31] ,  additional 
development  and  modifications  have  been  made  on  the  SAP  program 
as  follows: 

(1)  New  eigenvalue  solution  routines  using  a  determinant  search 
technique  combined  with  inverse  iteration  or  subspace 
iteration,  both  developed  by  Bathe  [26],  have  been  in- 
corporated into  the  program  for  the  determination  of 
frequencies  and  mode  shapes. 

(2)  In  the  dynamic  response  time  history  analysis  option 
of  the  program,  modifications  have  been  made  to  permit 
multiple  support  excitations  of  the  structural  system. 
One  may  wish  to  specify  this  type  of  input  for  some 
bridge  structures  which  are  extremely  long. 

(3)  A  subroutine  which  optimizes  the  block  storage  of  the 
stiffness  and  mass  matrices  of  the   system  on  the  basis 
of  available  core  storage  and  analysis  type  has  been 
incorporated  into  the  program. 

(4)  A  linear  elastic  circularly  curved  beam  element  has 

been  developed  and  added  to  the  element  library  for  model- 
ling curved  decks. 
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(5)  A  three-dimensional  boundary  spring  element  has  been 
added  to  the  element  library  for  modelling  foundation 
stiffnesses. 

(6)  A  linear  expansion  joint  element  has  been  developed  and 
added  to  the  element  library  for  modelling  joint 
behavior. 

Thus,  the  present  version  of  BSAP  can  be  used  to  analyze  linear 
bridge  structural  systems  modelled  by  any  combination  of  the 
following  element  types:   (1)  three-dimensional  truss  element, 
(2)  three-dimensional  straight  beam  element,  (3)  three-dimensional 
curved  beam  element,  (4)  three-dimensional  boundary  spring  element, 
and  (5)  linear  expansion  joint  element.   It  can  evaluate  (1)  response 
to  static  loadings,  (2)  frequencies  and  mode  shapes,  (3)  time  history 
responses  to  apply  dynamic  loadings  and/or  prescribed  ground 
excitations  (rigid  or  multiple),  and  (4)  modal  spectral  responses 
and  the  root-mean-square  values  using  prescribed  earthquake  spectra. 

The  organization  of  BSAP  follows  the  linear  analysis  pro- 
cedures established  previously  which  can  be  described  by  the  following 
operational  sequence  [31] : 

(1)  Nodal  coordinates  and  the  degrees  of  freedom  of  the 
system  are  generated  from  nodal  input  data.   The  nodal 
points  may  be  specified  as  fully  or  partially  con- 
strained in  any  of  its  displacement  coordinates.   Nodes 

at  multiple  support  input  points  are  identified  separately. 

(2)  Element  data  are  input;  then,  matrices  of  stiffness,  mass, 
and  local-to-global  coordinate  transformation  are 
generated  and  placed  in  low-speed  storage  units,  e.g., 
disc  files. 

(3)  Total  stiffness  and  mass  matrices  for  the  complete  system 
are  assembled  in  one  or  more  core  blocks  according  to 
element  connection  arrays .   Load  vectors  are  formed  in 
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the  case  of  static  analysis  and  the  support  coupling 
stiffness  matrix  £  is  assembled  in  the  case  of  time 
history  response  analysis  when  considering  multiple 
support  excitations. 

(4)  Equilibrium  equations  are  solved  for  static  response 

or  for  quasi-static  influence  matrix  B  when  prescribing 
multiple  support  excitations. 

(5)  For  dynamic  analysis,  a  specified  number  of  lowest 
frequencies  and  mode  shapes  are  determined. 

(6)  For  time  history  response  analysis,  applied  dynamic 
load  and/or  ground  excitation  time  functions  are  input 
and  digitized  at  equal  time  intervals.   The  response 
time  histories  are  calculated  using  mode  superposition 
procedure  as  described  previously. 

(7)  For  spectral  response  analysis,  earthquake  response  spectra 
are  input  and  modal  spectral  responses,  as  well  as  the 
root-mean-square  responses,  are  determined  based  on  pre- 
scribed modal  damping  ratios  and  the  computed  vibration 
frequencies. 

Computer  program  BSAP  is  coded  in  FORTRAN  IV.   All  storages 
are  allocated  in  terms  of  variable  dimensions  at  the  time  of 
execution,  and  an  out-core  block-by-block  solution  technique  is 
adopted.   Therefore,  the  program  is  capable  of  analyzing  large  bridge 
systems.   It  has  been  used  to  analyze  several  bridge  structures  using 
CDC  6600  and  7600  computers. 


B.   NONLINEAR  ANALYSIS  COMPUTER  PROGRAM  -  NEABS 

A  new  computer  program  NEABS  has  been  developed  for  earthquake 
analysis  of  nonlinear  bridge  structural  systems  of  the  type  considered 
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herein.   This  computer  program  uses  the  step-by-step  integration 
procedures  established  in  the  previous  chapter.   Either  the 
linear  acceleration  or  constant  acceleration  method  can  be  chosen 
for  integration.   The  excitations  can  result  from  applied  dynamic 
loads  and/or  from  support  motions  (rigid  or  multiple) .   The  program 
has  an  element  library  consisting  of  the  five  linear  elements  used 
in  BSAP  and  nonlinear  elements  as  follows:   (1)  three-dimensional 
elasto-plastic  flexural  column  (straight)  element  for  reinforced 
concrete  members,  (2)  three-dimensional  bilinear  boundary  spring 
element,  and  (3)  nonlinear  expansion  joint  element.   In  the  program, 
a  nonlinear  element  indicator  has  been  used  to  allow  a  system 
modelled  with  both  linear  and  nonlinear  elements. 

The  organization  of  computer  program  NEABS  can  be  described 
by  the  operational  sequence: 

(1)  Static  linear  equilibrium  equations  are  first  generated  and 
assembled  in  the  same  manner  as  in  the  linear  analysis 
program  BSAP.   They  are  then  solved  for  the  static  responses 
to  be  used  as  initial  conditions  in  the  subsequent  nonlinear 
earthquake  response  analysis.   In  this  same  process,  the 
initial  quasi-static  influence  matrix  is  also  determined, 

if  needed. 

(2)  Dynamic  loads  and/or  ground  excitation  time  functions  are 
input  and  are  digitized  for  equal  time  intervals. 

(3)  Checks  are  then  made  on  the  available  core  storage  to 
determine  whether  an  in-core  or  an  out-core  block-by-block 
solution  process  is  appropriate. 

(4)  For  an  in-core  solution,  an  in-core  step-by-step  solution 
package  is  used  to  determine  the  response  time  histories 
of  the  system  based  on  the  procedures  described  previously. 


-64- 


(5)   If  an  out-core  block-by-block  solution  is  necessary,  an 
out-core  step-by-step  solution  package  can  be  used  to 
evaluate  the  response  time  histories. 

Computer  program  NEABS  is  written  in  FORTRAN  IV  and  is  coded 
in  variable-dimension  form.   The  program  has  been  used  to  determine 
the  nonlinear  earthquake  responses  of  a  number  of  bridge  structures 
using  the  CDC  6600  and  7600  computers. 

Because  of  the  extensive  computational  effort  involved  in  a 
nonlinear  earthquake  response  analysis,  the  out-core  step-by-step 
solution  package  requires  many  slow-speed  storage  operations. 
Therefore,  this  procedure  can  be  very  costly  and  insufficient.. 
For  the  bridge  structures  analyzed  herein,  the  out-core  solution 
has  not  been  required. 
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VII   NUMERICAL  RESULTS 

The  previously  described  mathematical  modelling  and  analysis 
procedures  have  been  employed  to  determine  the  response  of  over  10 
bridge  structural  systems  using  computer  programs  BSAF  and  W.   jig. 
These  systems  are  basically  patterned  after  3  prototype  L. 
structures  but  with  variations  in  certain  p$fj§j^|  --  -      &&, 
3  prototype  structures  are  the  5/14  South  Connector  Overcrossing  of 
the  Golden  State  freeway  and  Antelope  Valley  freeway  interchange, 
the  curved  Figueroa  Street  Undercrossing  Connector  in  the  Los  Angeles 
area,  and  a  straight  version  of  the  Figueroa  Street  Undercrossing 
Connector.   The  structural  data  and  numerical  results  obtained  for  each 
bridge  system  are  presented  subsequently  along  with  summaries  of  extreme- 
values  of  response. 

A.   5/14  SOUTH  CONNECTOR  OVERCROSSING 

1.   The  Structural  System  -  The  structural  system  of  this  bridge 
consists  of  a  curved,  continuous,  reinforced  concrete  box  girder  deck 
of  9  spans  (Nos.  1-9  from  left  to  right  in  Fig.  14)  supported  on  8 
single,  central,  reinforced  concrete  columns  (Nos.  2-9  in  Fig.  14),  and 
2  reinforced  concrete  diaphragm  abutments  (Nos.  1  and  10  in  Fig.  14). 
The  deck  has  a  radius  of  curvature  of  667  ft.  measured  to  its  center 
line.   Four  expansion  joints  (1-4  from  left  to  right  in  Fig.  14)  divide 
the  deck  into  5  continuous  segments.  The  structural  layout  of  this 
bridge  system  is  shown  in  Fig.  14  and  the  structural  properties  of  the 
bridge  components  are  given  in  Table  1.  During  the  San  Fernando  earth- 
quake of  February  9,  1971,  this  bridge  suffered  severe  damage.  Central 
spans  Nos.  3  and  4  collapsed  along  with  column  No.  4. 

For  purpose  of  analysis,  this  bridge  has  been  idealized  by  a 
discrete-parameter  system  consisting  of  71  nodal  points  of  which  the 
10  located  at  the  bases  of  the  8  columns  and  the  2  abutments  are  fixed 
to  the  ground;  thus,  a  total  of  366  (61x6)  degrees  of  freedom  are 
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present  in  the  system  (Fig.  15) .  For  linear  analysis,  a  system  of 
28  three-dimensional  linear  straight  beam  elements,  38  three- 
dimensional  linear  curved  beam  elements,  and  4  linear  expansion 
joint  elements  have  been  used  to  model  the  structure  as  shown  in 
Fig.  16. 

For  nonlinear  analysis,  linear  elastic  straight  beam  elements 
Nos.  3-28  as  shown  in  Fig.  16  are  replaced  by  26  corresponding  three- 
dimensional  elasto-plastic  flexural  column  elements  and  the  4  linear 
expansion  joint  elements  are  replaced  by  nonlinear  elements. 

The  ultimate  column  strengths  as  computed  by  program  YIELD  are 
given  in  Table  2 .   The  corresponding  generalized  yield  function 
constants  are  shown  in  Table  3. 

For  nonlinear  analysis,  three  different  systems  of  longitudinal 
restrainers  have  been  considered  for  the  expansion  joints  (1)  a 
strong  system  consisting  of  3  short  tie  bars  each  having  a  2^  in. 
diameter,  a  4ft.  length,  and  a  yield  strength  of  480  kips,  (2)  a 
weak  system  consisting  of  3  short  tie  bars  each  having  a  1-  in. 
diameter,  a  4ft.  length,  and  a  yield  strength  of  70.8  kips,  and  (3)  a 
system  consisting  of  no  ties.   The  longitudinal  tie  bars  are  located 
across  each  expansion  joint  at  the  centers  of  the  3  box  girder  cells 
which  are  9ft.  apart.   Prior  to  the  San  Fernando  earthquake,  longitudinal 
restrainer  system  No.  2  was  generally  used.   However  following  this 
earthquake,  restrainer  system  No.  1  has  been  proposed  [32] . 

In  all  nonlinear  cases  studied,  the  vertical  restrainers  of  the 
expansion  joints  were  assumed  to  have  infinite  yield  strengths,  the 
Coulomb  coefficients  of  friction  were  set  at  0.4,  the  expansion  joint 

gaps  A  were  taken  equal  to  1  in.,  and  the  tie  gaps  were  assumed  equal 

G 

to  zero. 

Foundation  flexibilities  are  approximated  in  all  cases  by  assuming 
the  columns  and  abutments  to  be  fully  fixed  at  a  depth  10ft.  below  the 
ground  profile. 

2.  Mode  Shapes  and  Frequencies  -  Frequencies  and  corresponding 
mode  shapes  were  determined  for  two  structural  cases  using  program 
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BSAP.   In  the  first  case,  zero  friction  was  assumed  for  the  expansion 
.joints  and,  in  the  second  case,  infinite  friction  was  assumed.  The 
computed  frequencies  and  corresponding  periods  for  the  lowest  10 
modes  of  vibration  are  shown  in  Table  4  for  the  zero  friction  case  and 
in  Table  5  for  the  infinite  friction  case.   The  shapes  of  the  4 
lowest  modes  are  plotted  in  Figs.  17  and  18  for  the  zero  and  infinite 
friction  cases,  respectively. 

3.  Spectral  Responses  -  Spectral  responses  of  the  linear  bridge 
system  subjected  to  the  N-S  El  Centro  1940  earthquake  motions  in 
its  transverse  direction,  i.e.,  the  global  Y  direction  in  Fig.  15 
were  determined  for  both  the  zero  and  infinite  friction  cases  using 
program  BSAP.  A  damping  ratio  of  0.02  was  used  for  all  modes  in  these 
analyses. 

Spectral  responses  for  maximum  bending  moments  about  the  local  z 
axes  in  columns  Nos.  4,  5,  and  7  are  shown  for  the  lowest  10  modes 
in  tables  6  and  7  for  the  zero  and  infinite  friction  cases,  respectively 
Root-mean-square  values  for  all  10  modes  are  also  shown  in  these  tables. 

4.  Nonlinear  Earthquake  Responses  -  Nonlinear  earthquake  response 
analyses  using  computer  program  NEABS  were  carried  out  for  five 
structural  cases.   In  all  cases,  the  ground  excitation  inputs  were 
derived  from  an  artificial  earthquake  record  (OVH  accelerogram) 
generated  by  Chopra  [33]  as  shown  in  Fig.  19.   This  accelerogram  was 
generated  to  simulate  the  ground  motions  produced  by  the  San  Fernando 
earthquake  at  the  site  of  the  Olive  View  Hospital  located  about  6 
miles  southwest  of  the  epicenter.   It  has  a  peak  acceleration  of  0.5g, 

a  uniform  phase  of  high  intensity  shaking  for  8  seconds,  and  a  spectrum 
intensity  of  4.01  ft.  which  is  approximately  1.5  times  greater  than  the 
N-S  El  Centro  1940  earthquake. 

For  all  five  nonlinear  cases  studied,  rigid  ground  excitations 
were  assumed  and  the  intensity  levels,  i.e.,  peak  accelerations,  were 
adjusted  to  those  values  shown  in  Table  8.  The  structural  parameter 
varied  in  these  five  cases  was  a  type  of  longitudinal  restrainer  system. 
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The  numerical  results  obtained  for  these  five  cases  are  based  on- 
damping  constants  a  =  0.0419  and  3  =  0.0079  which  correspond  tc 
damping  ratios  of  0.02  in  the  first  two  transverse  modes  of  vibration 
for  the  zero  friction  linear  case  as  shown  in  Fig.  17.   In  all  cases, 
a  time  interval  of  0.02  seconds  and  the  constant  acceleration  method 
have  been  used  in  step-by-step  integration.   Representative  time 
histories  of  response  for  these  five  cases  are  presented  as  follows: 

Case  1  -  The  horizontal  accelerations  (absolute)  and  the 
displacements  in  the  global  X  and  Y  directions  at  the  top 
of  columns  Nos.  3,  4,  and  5  are  shown  in  Figs.  20,  21,  and 

22,  respectively.   The  vertical  accelerations  and  displace- 
ments at  the  center  of  spans  Nos.  3  and  4  are  given  in  Fig. 

23.  The  generalized  forces  at  the  center  of  span  No.  3  are 
shown  in  Fig.  24.   The  bending  moments  and  corresponding 
plastic  bending  rotations  at  the  bases  of  columns  Nos.  3, 

4,  5,  and  6  are  presented  in  Figs.  25,  26,  27,  and  28,  respective- 
ly.  Longitudinal  joint  separations  at  expansion  joints  Nos.  1  and 
2  are  shown  in  Fig.  29.   The  forces  and  corresponding  plastic 
elongations  for  the  3  longitudinal  tie  bars  at  expansion  joints 
Nos.  ]  and  2  are  given  in  Figs.  30  and  31,  respectively.   The 
forces  in  the  vertical  restrainers  of  expansion  joints  Nos.  1 
and  2  are  presented  in  Fig.  32. 

Case  2  -  The  horizontal  accelerations  and  displacements  at  the 
top  of  columns  Nos.  3,  4,  and  5  are  shown  in  Figs.  33,  34,  and 
35,  respectively.   The  vertical  accelerations  and  displacements 
at  the  centers  of  spans  Nos.  3  and  4  are  given  in  Fig.  36. 
The  generalized  forces  at  the  center  of  span  No.  4  are  shown 
in  Fig.  37.   The  bending  moments  and  corresponding  plastic 
bending  rotations  at  the  bases  of  columns  Nos.  3,  4,  5,  and  6 
are  presented  in  Figs.  38,  39,  40,  and  41,  respectively.   The 
longitudinal  joint  separations  at  expansion  joints  Nos.  1  and  2 
are  shown  in  Fig.  42.   The  longitudinal  tie  bar  forces  and 
corresponding  plastic  elongations  at  expansion  joints  Nos.  1 
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and  2  are  shown  in  Figs.  43  and  44,  respectively.  Vertical 
restrainer  forces  at  expansion  joints  Nos.  1  and  2  are  given 
in  Fig.  45. 

Case  3  -  This  particular  bridge  system  experienced  excessively 
large  amplitudes  of  response.  For  example,  note  the  large 
horizontal  displacements  which  occurred  at  the  top  of  column 
No.  4  as  shown  in  Fig.  46  and  the  large  joint  separations 
which  developed  at  expansion  joint  No.  2  as  presented  in  Fig.  47. 
Due  to  the  fact  that  very  large  geometry  changes  resulted  in 
this  case,  the  numerical  results  obtained  are  not  valid.   They 
do  however  strongly  indicate  that  some  portions  of  the  overall 
structure  would  collapse. 

Case  4  -  The  horizontal  accelerations  and  displacements  at 

the  top  of  column  No.  4  are  shown  in  Fig.  48.   The  vertical 

accelerations  and  displacements  at  the  top  of  this  same 

column  and  at  the  center  of  span  No.  4  are  presented  in  Fig.  49. 

The  generalized  forces  in  the  box  girder  at  the  center  of  span  No. 

4  are  given  in  Fig.  50.   The  bending  moments  and  corresponding 
plastic  bending  rotations  at  the  bases  of  columns  No.  3,  4,  and 

5  are  shown  in  Figs.  51,  52,  and  53,  respectively.   The  longi- 
tudinal joint  separations  at  expansion  joint  No.  2  are  presented 
in  Fig.  54.   The  longitudinal  tie  bar  forces  and  corresponding 
plastic  elongations  at  expansion  joint  No.  2  are  given  in  Fig.  55. 

Case  5  -  The  horizontal  accelerations  and  displacements  at  the 
top  of  column  No.  4  are  given  in  Fig.  56  while  the  vertical 
accelerations  and  displacements  at  the   top  of  this  same 
column  and  at  the  center  of  span  No.  4  are  presented  in  Fig.  57. 
The  generalized  forces  in  the  box  girder  at  the  center  of  span 
No.  4  are  shown  in  Fig.  58.   The  bending  moments  and  corresponding 
plastic  bending  rotations  at  the  bases  of  columns  Nos.  3,  4,  and 
5  are  presented  in  Figs.  59,  60,  and  61,  respectively.   The 
longitudinal  joint  separations  at  expansion  joint  No.  2  as 
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shown  in  Fig.  62  and  the  longitudinal  tie  bar  forces  and 
corresponding  plastic  elongations  are  shown  in  Fig.  63. 

B.   CURVED  FIGUEROA  STREET  UNDERCROSSING  CONNECTOR 

1.   The  Structural  System  -  The  structural  system  of  this 
bridge  consists  of  a  curved,  continuous,  reinforced  concrete  box 
girder  deck  of  6  spans  (Nos .  1-6  from  left  to  right  in  Fig.  64) 
supported  on  5  reinforced  concrete  columns  (Nos.  2-6  in  Fig.  64) 
and  2  reinforced  concrete  diaphragm  abutments  (Nos .  1  and  7  in 
Fig.  64).   The  deck  has  a  radius  of  curvature  of  886  ft.  measured 
to  its  center  line.   Only  one  expansion  joint  is  present  (No.  1). 
All  5  columns  have  a  length  of  40ft.  and  are  flared  at  their  upper 
ends  as  shown  in  Fig.  65.   These  columns  are  rigidly  connected  at 
their  bases  to  massive  reinforced  concrete  footings  (25'  x  25'  x  6') 
which  are  in  turn  supported  on  36  piles  each.  Therefore,  the 
foundation  stiffnesses  are  sufficiently  high  so  that  the  columns 
can  be  assumed  as  rigidly  fixed  at  their  bases.   The  box  girders  are 
rigidly  connected  to  the  diaphragm  abutments  which  are  relatively 
short  (6  ft.).  These  abutments  are  supported  on  spread  footings 
which  provide  high  shear  resistances  but  low  bending  resistances. 
The  structural  properties  of  this  bridge  system  are  summarized  in 
Table  9. 

The  discrete  parameter  system  used  to  model  this  bridge  is 
shown  in  Fig.  66.   It  contains  62  nodal  points  of  which  5  are 
located  at  the  fixed  bases  of  columns  and  2  are  located  at 
partially  fixed  bases  which  provide  full  constraint  in  the  three 
translational  directions  and  in  the  rotational  direction  about  a 
longitudinal  axis  parallel  to  the  deck;  however,  no  significant 
degree  of  fixity  is  provided  in  the  other  two  rotational  directions. 
Considering  these  support  constraints,  the  overall  system  has  274 
degrees  of  freedom. 

For  linear  analysis,  a  finite  element  system  consisting  of  27 
three-dimensional  straight  beam  elements,  23  three-dimensional 
curved  beam  elements,  and  1  linear  expansion  joint  element  as  shown 
in  Fig.  67  was  used.  Each  column  was  modelled  by  5  straight  beam 
elements  of  which  3  were  used  to  model  the  upper  flared  section 
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and  2  were  used  to  model  the  prismatic  portion.  For  nonlinear 

analysis,  the  2  straight  beam  elements  of  the  prismatic  section 

were  replaced  by  2  corresponding  elasto-plastic  flexural  column 

elements.  The  ultimate  strengths  of  these  elements  as  calculated 

by  YIELD  are  shown  in  Table  10.  The  corresponding  generalized  yield 

function  constants  are  given  in  Table  11. 

The  expansion  joint  of  this  bridge  uses  a  longitudinal  restrainer 

system  consisting  of  6  cable  units,  placed  across  the  joint  at 

distances  of  +  8.92,  +  5.85,  +  2.92  ft.  from  the  center  line  of  the 

deck.   Each  cable  unit  consists  of  5  -  lr  in.  diameter  cables,  127.5  ft. 

long,  having  yield  strengths  of  154.6  kips  each.  The  expansion  joint 

is  constructed  with  an  initial  joint  gap  A  =  6  in.  but  with  all  tie 

G 

gaps  A  equal  to  zero.  The  vertical  restrainers  at  the  expansion 
joint  are  assumed  to  have  infinite  yield  strengths.   The  Coulomb 
friction  coefficient  is  again  assumed  to  be  0.4.   The  expansion  joint 
width  d  is  taken  as  27ft. 

The  longitudinal  restrainer  system  used  in  this  bridge  is  a  new 
version  recently  developed  by  the  California  State  Division  of 
Highways  [32] . 

2.  Mode  Shapes  and  Frequencies  -  The  mode  shapes  and  frequencies  of 
this  bridge  system  have  been  determined  only  for  the  case  of  zero 
friction  in  the  expansion  joint.  The  computed  frequencies  and  correspond- 
ing periods  for  the  lowest  10  modes  are  listed  in  Table  12.  The  shapes 
of  the  lowest  4  modes  are  plotted  in  Fig.  68. 

3.  Nonlinear  Earthquake  Responses  -  The  nonlinear  response  of  this 
bridge  system  has  been  determined  using  NEABS  for  three  different  cases. 
In  all  three  cases,  rigid  ground  excitations  corresponding  to  the  OVH 
accelerogram  were  used  as  shown  in  Fig.  19.  The  structural  parameter 
varied  in  these  studies  were  the  yield  strengths  of  the  columns.   Both 
elastic  and  elasto-plastic  flexural  columns  were  considered.   Table  13 
describes  the  column  types  used  in  each  case  and  also  shows  the 
intensities  of  excitations  which  were  prescribed. 

Numerical  results  have  been  obtained  for  all  three  cases  using 


-72- 


damping  constants  a  =  0.09  and  3  =  0.0012  which  correspond  to  a 
damping  ratio  of  0.02  for  the  lowest  2  modes  of  vibration  for  the 
zero  friction  case.   In  all  3  cases,  a  time  interval  of  0.02  seconds 
has  been  used  in  step-by-step  integration.   Selected  time  history 
responses  are  presented  as  follows: 

Case  1  -  The  horizontal  accelerations  and  displacements 
at  the  top  of  column  No.  3  and  4  are  presented  in  Figs.  69 
and  70,  respectively.   The  vertical  accelerations  and  dis- 
placements at  the  top  of  column  No.  3  and  at  the  center  of 
span  No.  3  are  shown  in  Fig.  71.   The  generalized  forces  in 
the  box  girder  at  the  center  of  span  No.  3  are  given  in 
Fig.  72  while  the  elastic  bending  moments  at  the  bases  of 
columns  Nos.  3  and  4  are  presented  in  Fig.  73.   The 
longitudinal  joint  separations  are  shown  in  Fig.  74  and 
the  longitudinal  cable  forces  at  the  joint  are  given  in  Fig.  75. 
These  cable  forces  remain  elastic  during  the  entire  response 
history. 

Case  2  -  The  horizontal  accelerations  and  displacements  at 
the  top  of  columns  Nos.  3  and  4  are  shown  in  Figs.  76  and  77, 
respectively.  The  vertical  accelerations  and  displacements 
at  the  top  of  column  No.  3  and  at  the  center  of  span  No.  3 
are  given  in  Fig.  78.   The  generalized  forces  in  the  box  girder 
at  the  center  of  span  No.  3  are  presented  in  Fig.  79  while  the 
bending  moments  and  corresponding  plastic  rotations  at  the 
bases  of  columns  Nos.  3  and  4  are  shown  in  Figs.  80  and  81, 
respectively.  The  longitudinal  joint  separations  are  given 
in  Fig.  82  and  the  longitudinal  cable  forces  are  presented  in 
Fig.  83.  These  cable  forces  also  remain  elastic  during  the 
entire  response  history. 

Case  3  -  The  horizontal  accelerations  and  displacements  at 
the  top  of  column  No.  4  are  shown  in  Fig.  84  while  the  vertical 
accelerations  and  displacements  at  the  top  of  column  No.  3  and 
at  the  center  of  span  No.  3  are  given  in  Fig.  85.  The 
generalized  forces  in  the  box  girder  at  the  center  of  span  No.  3 
are  presented  in  Fig.  86  while  the  bending  moments  and  correspond - 
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ing  plastic  rotations  at  the  bases  of  columns  Nos.  3  and  4 
are  shown  in  Figs.  87  and  88,  respectively.   The  longitudinal 
joint  separations  are  presented  in  Fig.  89  and  the  longi- 
tudinal cable  forces  are  given  in  Fig.  90.   It  is  significant 
to  note  that  the  cables  remained  elastic  throughout  the  time 
history  of  response.       : 

C.   STRAIGHT  FIGUEROA  STREET  UNDERCROSSING  CONNECTOR 

1.  The  Structural  System  -  The  structural  system  of  this  bridge 
is  identical  to  that  of  the  curved  Figueroa  Street  Undercrossing 
Connector  described  previously  except  that  the  curvature  of  the 
deck  has  been  removed.    The  general  layout  of  the  bridge  system  is 
shown  in  Fig.  91  and  the  discrete  parameter  system  selected  for  analysis 
is  given  in  Fig.  92.   This  discrete  parameter  system  has  the  same 
nodal  point  arrangement  as  that  selected  for  the  curved  bridge  model. 
For  linear  analysis,  the  structure  has  been  modelled  by  50  three- 
dimensional  straight  beam  elements  and  1  linear  expansion  joint  as 
shown  in  Fig.  93.   For  nonlinear  analysis,  2  three-dimensional  elasto- 
plastic  flexural  elements  have  been  used  to  model  the  lower  22ft. 
prismatic  section  of  each  column.   The  ultimate  strengths  and  cor- 
responding yield  function  constants  for  these  columns  are  listed  in 
Tables  10  and  11. 

2.  Mode  Shapes  and  Frequencies  -  Mode  shapes  and  frequencies  of 
this  bridge  system  have  been  determined  for  the  zero  friction  case 
only.   Frequencies  and  corresponding  periods  for  the  lowest  10 
modes  are  listed  in  Table  14  while  the  corresponding  shapes  of  the 
lowest  4  modes  are  shown  in  Fig.  94. 

3.  Nonlinear  Earthquake  Responses  -  Nonlinear  response  analyses 
for  this  system  have  been  carried  out  using  NEABS  for  only  one  case. 
The  rigid  ground  excitations  used  corresponded  with  the  OVH  accelero- 
grams.  The  intensities  used  were  0.5g  in  the  transverse  direction 
(global  Y  direction)  and  0.3g  in  the  vertical  direction  (global  Z 
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direction) . 

The  nonlinear  analysis  was  carried  out  using  damping  constants 
a  =  0.09  and  $  =  0.0013  which  correspond  to  a  damping  ratio  of 
0.02  for  the  lowest  2  modes  of  vibration  in  the  zero  friction  case, 
A  time  interval  of  0.02  seconds  has  also  been  used.   The  horizontal 
accelerations  and  displacements  at  the  top  of  columns  Nos.  3  and  4 
are  shown  in  Figs.  95  and  96,  respectively,  while  the  vertical 
accelerations  and  displacements  at  the  top  of  column  No.  3  and  at 
the  center  of  span  No.  3  are  given  in  Fig.  97.   The  generalized 
forces  in  the  box  girder  at  the  center  of  span  No.  3  are  presented 
in  Fig.  98  and  the  bending  moments  and  corresponding  plastic  rotations 
at  the  base  of  columns  Nos.  3  and  4  are  shown  in  Figs.  99  and  100, 
respectively.   The  longitudinal  joint  separations  are  given  in 
Fig.  101  and  the  longitudinal  cable  forces  are  shown  in  Fig.  102. 
It  should  be  noted  that  the  cable  forces  remained  elastic  through- 
out the  time  history  of  response. 

D.   SUMMARY  OF  EXTREME -VALUES  OF  NONLINEAR  RESPONSES 

The  extreme -values  of  the  responses  previously  presented  for  the 
5  cases  of  the  5/14  South  Connector  Overcrossing,  the  3  cases  of  the 
curved  Figueroa  Street  Undercrossing  Connector,  and  the  single  case 
of  the  straight  Figueroa  Street  Undercrossing  Connector  are  summarized 
as  follows: 

1 .  Maximum  Accelerations  and  Displacements  -  The  maximum 
horizontal  accelerations  (absolute)  and  displacements  at  the  top 
of  columns  No.  2-7  are  summarized  in  Tables  15  and  16,  respectively, 
for  the  5/14  South  Connector  Overcrossing.   The  maximum  values 
shown  are  the  vectorial  sum  of  the  horizontal  components  measured 
along  the  global  X  and  Y  axes.   Likewise,  the  maximum  horizontal 
accelerations  and  displacements  at  the  top  of  each  column  for  the 
3  cases  of  the  curved  and  the  single  case  of  the  straight  Figueroa 
Street  Undercrossing  Connector  are  shown  in  Tables  17  and  18, 
respectively. 
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The  maximum  vertical  accelerations  and  displacements  at  the 
center  of  spans  are  shown  in  Tables  19  and  20,  respectively,  for 
the  5/14  South  Connector  Overcrossing  and  are  shown  in  Tables 
21  and  22,  respectively,  for  the  curved  and  straight  Figueroa 
Street  Undercrossing  Connector. 

2.  Locations  of  Maximum  Flexural  Yielding  in  Columns  -  The 
locations  of  maximum  flexural  yielding  in  the  columns  of  the  5/14 
South  Connector  Overcrossing  are  shown  in  Fig.  103  for  all  5  cases 
investigated.   Similar  locations  for  Cases  2  and  3  of  curved 
Figueroa  Street  Undercrossing  Connector  and  the  single  case  of  the 
straight  Figueroa  Street  Undercrossing  Connector  are  shown  in 
Fig.  104.   It  should  be  noted  that  the  hinge  locations  in  these 
figures  represent  locations  of  maximum  flexural  yielding  occurring 
during  the  first  10  seconds  of  response. 

3.  Maximum  Local  Bending  Ductility  Factors  at  the  Bases  of 
Columns  -  The  maximum  local  bending  ductility  factor  U  at  a 
yield  hinge  of  a  column  is  defined  as  the  ratio  of  the  maximum 
flexural  rotation  to  the  yield  rotation  0  ,  i.e. 


ep  +  eY 


(113) 


where  6  is  the  plastic  component  of  flexural  rotation.  For  a 

column  of  finite  dimension,  yielding  will  not  occur  in  concentrated 

hinges;  rather,  it  will  occur  over  finite  lengths  h  which  are 

approximately  equal  to  the  transverse  dimensions  of  the  column 

cross-section  in  the  directions  normal  to  the  axes  of  rotation 

y 
[20] .   Therefore,  the  flexural  yield  rotation  0  can  be  estimated 


by  the  relation 


M 
0Y  =  -£  h  (114) 
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where  M  is  the  ultimate  bending  moment  which  in  the  present 

analysis  is  assumed  equal  to  the  yield  bending  moment. 

Y       Y 
The  flexural  yield  rotations  6   and  6   about  the  local  y  and 
2  y       z 

z  axes  of  the  columns  of  the  5/14  South  Connector  Overcrossing  have 
been  computed  using  Eq.  (114)  and  the  results  are  listed  in  Table  23. 
Similar  yield  rotations  for  the  columns  of  the  Figueroa  Street 
Undercrossing  Connector  are  shown  in  Table  24. 

The  maximum  local  bending  ductility  factors  at  the  bases 
of  the  columns  of  the  5/14  South  Connector  Overcrossing  are  computed 
using  Eq.  (113)  for  Cases  1,  2,  4,  and  5.  The  results  are  shown 
in  Table  25.   Similar  ductility  factors  for  Cases  2  and  3  of  the 
curved  Figueroa  Street  Undercrossing  Connector  and  the  single  case 
of  the  straight  Figueroa  Street  Undercrossing  Connector  are  given 
in  Table  26.  The  values  presented  in  Tables  25  and  26  are  absolute 
maxima  computed  for  the  bending  rotations  about  the  local  y  and  z  axes. 

4.  Maximum  Longitudinal  Joint  Separations  and  Longitudinal 
Re strainer  Tie  Ductility  Factors  -  The  maximum  longitudinal  expansion 
joint  separations  at  the  center  lines  of  the  bridge  decks  are  shown 
in  Table  27  for  expansion  joints  Nos.  1-4  in  Cases  1,  2,  4,  and  5 
of  the  5/14  South  Connector  Overcrossing.  Similar  joint  separations 
are  shown  in  Table  28  for  the  single  expansion  joint  of  the  Figueroa 
Street  Undercrossing  Connector  (3  cases  for  the  curved  undercrossing  and 
1  case  for  the  straight  crossing)  . 

The  longitudinal  restrainer  tie  ductility  factor  yT  is  defined  as 
the  ratio  of  the  maximum  tie  elongation  to  its  yield  elongation,  i.e. 

E    P 
u   +  u 

T      T 

UT-  -i— -i-  (115) 

T 

where  uT?  and  uTE  are  defined  by  Eq.  (53)  of  Chapter  IV.  Faster  UT 
is  computed  in  each  case  for  all  longitudinal  restrainer  ties  at  the 
expansion  joints.  The  maximum  vslues  of  y_  for  eaeh  expaniien  joint 
are  listed  in  Table  27  for  4  cases  of  the  5/14  South  Cenneefcer 
Overcrowing  and  in  Table  28  for  3  cases  of  the  curved  Figueroa  Street 
Underleasing  Connector  and  for  the  single  ease  of  the  straight 
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Figueroa  Street  Undercrossing  Connector, 
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Table  2    Ultimate  Strength  of  the  Columns  of 
5/14  South  Connector  Overcrossing 


Column 
No. 

(  psi  ) 

Reinforce- 
ment 

#18  bars 

Po 
10  *  (k) 

Myo 
10 k    (k-ft) 

Mzo 
10  **  (k-ft) 

#  2 

3,500 

40 

,2.579 

1.386 

2.537 

#  3,  5 

3,500 

40 

2.948 

1.671 

2.598 

#  4 

4,000 

52 

3.470 

2.156 

3.331 

#  6 

4,000 

40 

2.465 

1.102 

2.517 

#  7,  8 

4,500 

48 

2.821 

1.320 

2.997 

#  9 

3,500 

48 

2.365 

1.303 

2.900 

Table  3    Generalized  Yield  Function  Constants  of  the  Columns  of 
5/14  South  Connector  Overcrossing 


Column 
No. 

a 

ai 

a2 

aa 

b 

bi 

b2 

b3 

#  2 

2.0 

-3.10 

-3.83 

0.273 

2.0 

-2.97 

-4.21 

-0.244 

#  3 

2.0 

-3.63 

-4.53 

0.102 

2.0 

-3.55 

-4.80 

-0.245 

#  4 

2.0 

-3.20 

-4.09 

0.101 

2.0 

-3.12 

-4.34 

-0.225 

#  6 

2.0 

-2.98 

-3.44 

0.546 

2.0 

-2.78 

-4.00 

-0.219 

#  7,  8 

2.0 

-2.81 

-3.32 

0.489 

2.0 

-2.61 

-3.85 

-0.238 

#  9 

2.0 

-2.22 

-2.79 

0.422 

2.0 

-2.01 

-3.24 

-0.232 
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Table  4    Frequencies  and  Periods  of  5/14  South 
Connector  Overcrossing  -  Zero  Friction 


Mode 
No. 

Frequency 
(rad  /  sec) 

Period 
(sec) 

1 

1.226 

5.124 

2 

2.243 

2.801 

3 

3.121 

2.013 

4 

5.549 

1.132 

5 

7.929 

0.924 

6 

9.031 

0.696 

7 

9.123 

0.689 

8 

9.507 

0.661 

9 

12.241 

0.513 

10 

12.902 

0.487 

Table  5    Frequencies  and  Period  of  5/14  South 

Connector  Overcrossing  -  Infinite  Friction 


Mode 
No. 

Frequency 
(rad  /  sec) 

Period 
(sec) 

1 

4.500 

1.391 

2 

6.209 

1.012 

3 

9.109 

0.690 

4 

9.585 

0.655 

5 

11.301 

0.556 

6 

13.487 

0.466 

7 

14.262 

0.441 

8 

14.668 

0.428 

9 

15.123 

0.415 

10 

16.013 

0.392 
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Table  6    Spectral  Responses  of  5/14  South  Connector 

Overcrossing  to  the  N-S  El  Centro  1940  Earthquake 
— Zero  Friction  Case 


Mode 
No. 

Maximum  Bending  Moments  in  Columns  * 
About  the  Local  z  Axis  (k-ft) 

Column  #  4 

Column  #  5 

Column  #  7 

1 

22.5 

102.7 

1.5 

2 

29,311.0 

31,334.0 

2,115.4 

3 

2.1 

398.7 

30.4 

4 

1,635.1 

7,828.6 

2,455.0 

5 

639.7 

2,359.4 

17,458.0 

6 

902.3 

5,519.3 

12,080.0 

7 

40.0 

18.5 

16.9 

8 

747.9 

762.2 

82.9 

9 

437.8 

697.4 

1,664.9 

10 

1,841.6 

791.3 

94.2 

R.M.S. 

29,448.0 

32,879.0 

21,540.0 

*  Maximum  moments  shown  all  occur  at  the  base  of  each  column 
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Table  7    Spectral  Responses  of  5/14  South  Connector  Overcrossing 
to  the  N-S  El  Centre  1940  Earthquake  —  Infinite 
Friction  Case 


Mode 
No. 

Maximum  Bending  Moments  in  Columns  ** 
about  the  Local  z  Axis  (k-ft) 

Column  #  4 

Column  #  5 

Column  #  7 

1 

554.2 

480.7  * 

884.6 

2 

11,573.0 

10,141.0 

1,928.9  * 

3 

41.8  * 

446.9  * 

237.5 

4 

2,136.2 

4,358.1  * 

12,570.0 

5 

5.2 

5.8 

4.5 

6 

23.6 

95.5 

20.9  * 

7 

10.0 

1.4 

164.9  * 

8 

237.1 

12.8 

745.4 

9 

100.5 

79.8  * 

293.3 

10 

61.8 

38.9 

222.5 

R.M.S. 

11,785.0 

i            . 

10,574.0 

12,684.0 

*  Maximum  moment  occurs  at  the  top  of  the  column. 
**  All  valves  shown  are  maxima  at  the  base  of  the  columns,  except 
those  marked  by  *. 
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Table  10    Ultimate  Strength  of  the  Columns  of  Figueroa 
Street  Under crossing  Connector 


Column 
No. 

f  ' 
^c 

(lb/in2) 

Long. 
Bars 
#11 

*0 
10"  (k) 

10*  (k-ft) 

Mzo 
10**  (k-ft) 

#  2 

3,250 

96 

2.698 

2.456 

2.456 

#  3,4 

3,250 

87 

2.614 

2.253 

2.253 

#  5 

3,250 

75 

2.502 

1.978 

1.978 

#  6 

3,250 

100 

3.189 

3.555 

3.555 

Table  11    Yield  Function  Constants  of  the  Columns  of  Figueroa 
Street  Under crossing  Connector 


Column 
No. 

a 

ai 

a2 

a3 

b 

*>1 

b2 

*>3 

#  2 

2.0 

-1.86 

-3.29 

-0.440 

2.0 

-1.86 

-3.29 

-0.440 

#  3,4 

2.0 

-1.82 

-2.47 

0.352 

2.0 

-1.82 

-2.47 

0.352 

#  5 

2.0 

-2.44 

-3.89 

-0.455 

2.0 

-2.44 

-3.89 

-0.455 

#  6 

2.0 

-1.16 

-2.47 

-0.317 

2.0 

-1.16 

-2.47 

-0.317 
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Table  12    Frequencies  and  Periods  of  Figueroa 
Street  Undercrossing  Connector  — 
Zero  Friction  at  Expansion  Joint 


Mode 
No. 

Frequency 
(rad/sec) 

Period 
(sec) 

1 

7.072 

0.889 

2 

9.208 

0.682 

3 

10.075 

0.624 

4 

11.866 

0.530 

5 

13.924 

0.451 

6 

15.382 

0.409 

7 

20.080 

0.313 

8 

22.986 

0.273 

9 

24.171 

0.260 

10 

26.921 

0.233 
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Table  14    Frequencies  and  Periods  of  Straight 

Figueroa  Street  Undercrossing  Connector 
—  Zero  Friction  at  Expansion  Joint 


Mode 
No. 

Frequency 
(rad/sec) 

Period 
(sec) 

1 

7.213 

0.871 

2 

9.315 

0.675 

3 

10.098 

0.622 

4 

13.921 

0.451 

5 

13.965 

0.450 

6 

15.367 

0.409 

7 

20.315 

0.309 

8 

23.390 

0.269 

9 

26.481 

0.237 

10 

28.081 

0.224 
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Table  15    Maximum  Horizontal  Accelerations  (Absolute)  at  the  Top 
of  Columns  of  5/14  South  Connector  Overcrossing 


^^Column 
Case  NoN. 

#2 
(g) 

#3 

(g) 

(g) 

#5 

(g) 

#6 

(g) 

#7 

(g) 

Case  1 

0.93 

0.86 

0.62 

0.87 

0.83 

0.64 

Case  2 

0.94 

1.96 

1.01 

0.99 

1.81 

0.85 

Case  4 

0.44 

0.64 

0.36 

0.88 

0.53 

0.40 

Case  5 

1.01 

0.83 

0.70 

0.85 

0.72 

0.51 

Table  16    Maximum  Horizontal  Displacements  at  the  Top 

of  Columns  of  5/14  South  Connector  Overcrossing 


N.  Col. 
Case  NoS^ 

— — - - — — — ■ 

#2 
(ft) 

#3 

(ft) 

#4 

(ft) 

#5 

(ft) 

#6 
(ft) 

#7 
(ft) 

Case  1 

1.21 

2.50 

2.57 

1.96 

1.91 

1.03 

Case  2 

1.45 

3.32 

3.70 

3.07 

1.30 

0.90 

Case  4 

0.69 

1.41 

1.69 

1.47 

0.68 

0.48 

Case  5 

0.78 

1.69 

1.94 

1.80 

1.28 

0.74 
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Table  17    Maximum  Horizontal  Accelerations  (Absolute)  at  the  Top 
of  Columns  of  Figueroa  Street  Undercrossing  Connector 


N^Column 
Case  NoV. 

#2 

(g) 

#3 

(g) 

#4 

(g) 

#5 

(g) 

#6 

(g) 

Case  1 

1.01 

0.94 

1.38 

1.06 

0.68 

Case   2 

0.75 

0.74 

0.85 

0.76 

0.58 

Case   3 

0.49 

0.51 

0.55 

0.39 

0.44 

Straight 

0.86 

0.84 

0.91 

0.61 

0.72 

Table  18    Maximum  Horizontal  Displacements  at  the  Top  of 

Columns  of  Figueroa  Street  Undercrossing  Connector 


^\Column 
Case  NoV^ 

#2 
(ft) 

#3 
(ft) 

#4 
(ft) 

#5 
(ft) 

#6 
(ft) 

Case  1 

0.17 

0.32 

0.68 

0.61 

0.31 

Case   2 

0.13 

0.23 

0.56 

0.40 

0.23 

Case   3 

0.11 

0.15 

0.33 

0.26 

0.15 

Straight 

0.18 

0.34 

0.91 

0.68 

0.28 
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Table  19    Maximum  Vertical  Accelerations  (Absolute)  at  the  Center 
of  Spans  of  5/14  South  Connector  Overcrossing 


n.  Span 
Case  NoV^ 

#2 

(g) 

#3 

(g) 

#4 

(g) 

#6 

(g) 

#7 

(g) 

#8 

(g) 

Case  1 

0.71 

0.59 

0.62 

0.57 

0.73 

0.53 

Case  2 

0.71 

0.62 

0.59 

0.54 

0.71 

0.44 

Case  4 

0.41 

0.38 

0.29 

0.30 

0.44 

0.28 

Case  5 

0.21 

_ 

0.15 

0.16 

0.25 

0.16 

0.12 

Table  20    Maximum  Vertical  Displacements  at  the  Center 
of  Spans  of  5/14  South  Connector  Overcrossing 


^\Span 
Case  No>\ 

#2 
(ft) 

#3 
(ft) 

#4 
(ft) 

#6 
(ft) 

#7 
(ft) 

#8 
(ft) 

Case  1 

0.091 

0.106 

0.138 

0.047 

0.051 

0.084 

Case  2 

0.071 

0.118 

0.135 

0.096 

0.073 

0.085 

Case  4 

0.055 

0.076 

0.085 

0.047 

0.034 

0.059 

Case  5 
( 

0.055 

0.047 

0.054 

0.042 

0.022 

0.029 
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Table  21    Maximum  Vertical  Accelerations  (Absolute)  at  the  Center 
of  Spans  of  Figueroa  Street  Undercrossing  Connector 


^s.    Span 
Case  No>\ 

#2 

(g) 

#3 

(g) 

#4 

(g) 

#5 

(g) 

#6 

(g) 

Case  1 

0.90 

0.77 

1.60 

1.70 

0.95 

Case   2 

1.00 

0.58 

1.65 

1.65 

0.86 

Case   3 

0.55 

0.39 

0.97 

0.99 

0.56 

Straight 

1.12 

0.92 

1.30 

2.64 

0.99 

Table  22    Maximum  Vertical  Displacement  at  the  Center  of 

Spans  of  Figueroa  Street  Undercrossing  Connector 


^\Span 
Case  No>vs^ 

#2 
(ft) 

#3 
(ft) 

#4 
(ft) 

#5 
(ft) 

#6 
(ft) 

Case  1 

0.046 

0.129 

0.117 

0.262 

0.253 

Case  2 

0.051 

0.123 

0.161 

0.254 

0.238 

Case  3 

0.026 

0.077 

0.093 

0.147 

0.152 

Straight 

0.054 

0.155 

0.153 

0.375 

0.241 

• 
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Table  25    Maximum  Local  Bending  Ductility  Factors  at  the  Base 
of  Columns  of  5/14  South  Connector  Overcrossing 


^\Base  of 
Case  No/\^ 

Column 
#  2 

Column 
#  3 

Column 
#  4 

Column 
#  5 

Column 
#  6 

Column 

#  7 

Case  1 

1.3 

2.7 

1.3 

3.2 

4.6 

3.6 

Case  2 

4.8 

4.0 

3.0 

7.4 

12.2 

2.2 

Case  4 

1.5 

1.2 

1.0 

3.3 

2.5 

1.8 

Case  5 

1.6 

1.8 

1.1 

3.3 

7.8 

2.3 

Table  26    Maximum  Local  Bending  Ductility  Factors  at  the  Base 

of  Columns  of  Figueroa  Street  Undercrossing  Connector 


^\Base  of 
Case  No?Sss. 

Column 
#  2 

Column 
#  3 

Column 
#  4 

Column 
#  5 

Column 
#  6 

Case  1 

_* 

- 

- 

- 

- 

Case  2 

1.3 

4.0 

5.4 

4.7 

- 

Case  3 

- 

2.4 

4.5 

3.1 

- 

Straight 

2.3 

4.2 

22.2 

18.0 

1.2 

*  No  yielding  occurs. 
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Table  27  Maximum  Joint  Separations  and  Longitudinal  Restrainer 
Tie  Ductility  Factors  at  the  Expansion  Joints  of  5/14 
South  Connector  Overcrossing 


Case  No 

Max.  Separation  (ft) 

I 

Tie  Ductility  Factor 

EJ  #1 

EJ  #2 

EJ  #3 

EJ  #4 

EJ  #1 

EJ  #2 

EJ  #3 

EJ  #4 

Case  1 

i 

1.21 

0.92 

0.40 

0.21 

78.8 

58.4 

26.2 

13.4 

Case  2 

1.40 

1.63 

0.37 

0.77 

278.0 

33S.0 

94.5 

145.1 

Case  4 

0.79 

0.48 

0.07 

0.34 

134.0 

101.0 

28.3 

63.6 

Case  5] 

0.94 

0.63 

0.40 

0.22 

60.4 

41.2 

25.4 

14.4 

Table  28    Maximum  Joint  Separations  and  Longitudinal 
Restrainer  Tie  Ductility  Factors  at  the 
Expansion  Joint  of  Figueroa  Street  Undercrossing 
Connector 


Case  No. 

Joint  Separation 
(ft) 

Tie  Ductility 
Factor 

Case  1 

0.305 

_* 

Case  2 

0.237 

- 

Case  3 

0.133 

- 

Straight 

0.001 

*  No  yielding  occurs. 
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Fig.  24  Generalized  Forces  in  the  Girder  at  the  Center 
of  Span  #  4 
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Fig.  103  Location  of  Maximum  Flexural  Yielding  in  Columns  of 
5/14  South  Connector  Overcrossing 
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VIII   DISCUSSION  OF  NUMERICAL  RESULTS 


A.   SEISMIC  RESPONE  CHARACTERISTICS  OF  BRIDGE  SYSTEMS 

The  analytical  results  previously  presented  for  the  5/14 
South  Connector  Overcrossing  showed  that  this  bridge  experienced 
large  amplitude  response  when  subjected  to  severe  ground  shaking 
(Cases  1,  2,  3,  and  5).   For  Case  2  which  uses  weak  short  tie 
bars ,  the  longitudinal  deck  separations  exceeded  the  supporting 
ledge  dimension  (1.25  ft.)  at  expansion  joints  Nos.  1  and  2;  see 
Table  27.   The  magnitudes  of  these  separations  correspond  to  tie 
bar  ductility  factors  far  in  excess  of  those  actually  provided?  thus, 
early  failures  of  these  bars  would  result  allowing  high  amplitude 
response  similar  to  Case  3  (no  tie  bars)  to  develop.  This  high  ampli- 
tude response  would  result  in  collapse  of  the  structure. 

For  Cases  1  and  5  which  use  strong  short  tie  bars,  the  joint 
separations  never  exceeded  the  supporting  ledge  dimension  even  under 
conditions  of  severe  shaking;  thus,  collapse  of  the  structure  in 
these  2  cases  would  not  occur  if  the  required  tie  ductilities  are 
provided. 

The  analytical  results  presented  for  the  curved  Figueroa  Street 
Undercrossing  Connector  (Case  2)  showed  relatively  high  amplitude 
acceleration  response  for  the  condition  of  severe  excitations. 
Although  this  response  causes  some  f lexural  damages  at  the  bases  of 
columns  supporting  the  central  spans ,  the  oscillatory  motions  were 
quite  stable.  Due  to  the  fact  that  the  columns  are  relatively  short 
and  only  one  expansion  joint  is  present  in  the  deck,  the  deck  dis- 
placements and  the  expansion  joint  separations  are  relatively  small. 
The  longitudinal  restrainer  cables  at  the  joint  never  reach  their 
yield  condition  (0.64  ft.  elongation)  even  under  severe  shaking. 
The  maximum  joint  separation  in  this  case  was  0.3  ft. 

The  horizontal  time  histories  at  the  tops  of  all  columns  for 
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all  bridge  structures  examined  show  that  all  segments  of  deck  between 
expansion  joints  move  essentially  as  rigid  bodies  in  the  horizontal 
directions,  i*e,  transverse  bending  effects  in  the  deck  are  small. 
Therefore,  the  deck  segments  vibrate  horizontally  almost  as  rigid 
bodies  linked  together  at  the  expansion  joints.   The  principal  mode 
of  vibration  resulting  from  transverse  excitation  is  similar  to  the 
transverse  fundamental  mode  for  the  linear  zero  friction  case,  i.e., 
similar  to  mode  No.  2  shown  in  Fig.  17  ior  the  5/14  South  Connector 
Overcrossing  and  mode  No.  1  shown  in  Figs.  68  and  94  for  the  curved 
and  straight  Figueroa  Street  Undercrossing  Connector,  respectively. 

It  should  be  noted  that  for  curved  bridges,  the  transverse 
fundamental  mode  of  vibration  caused  large  forces  to  develop  in  the 
restrainer  tie  bars.   Therefore,  the  properties  of  the  restrainer  tie 
bars  greatly  affect  the  transverse  seismic  responses  of  such  bridges. 
These  effects  can  be  observed  by  comparing  the  results  for  Case  1   and 
2  of  the  5/14  South  Connector  Overcrossing. 

The  overall  horizontal  displacement  responses  for  all  cases  of 
the  5/14  South  Connector  Overcrossing  indicate  drifts  of  the  displace- 
ments in  directions  away  from  the  centers  of  curvature  of  the  decks, 
i.e.,  outward  motions.   This  characteristic  behavior  is  due  to  (1)  the 
relatively  high  stiffnesses  which  develop  as  expansion  joints  close 
with  inward  displacements  causing  the  deck  to  act  as  a  continuous  stiff 
arch  from  abutment  to  abutment,  and  (2)  the  relatively  low  stiffnesses 
which  remain  with  outward  displacements.   For  outward  displacements, 
the  stiffnesses  are  contributed  only  by  the  actions  of  columns  and 
restrainer  tie  bars. 

The  drifts  of  horizontal  displacements  in  the  outward  directions 
are  particularly  noticeable  for  those  cases  of  loading  which  cause 
plastic  deformations  to  occur  in  the  columns  and  in  the  restrainer  tie 
bars.   Because  of  this  biased  behavior,  the  plastic  deformations  tend 
to  accumulate  in  the  outward  direction.   As  a  result,  the  entire  bridge 
structure  acts  as  a  "shake-down"  system. 

In  all  cases  of  the  curved  Figueroa  Street  Undercrossing  Connector, 
the  overall  horizontal  displacement  response  was  oscillatory  and 
approximate3.y  equal  in  both  inward  and  outward  directions.   The  plastic 
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flexural  rotation  at  the  column  bases  was  also  oscillatory  and  did 
not  accumulate  in  the  outward  direction  as  in  the  case  of  5/14  South 
Connector  Overcrossing.   The  stable  response  of  the  curved  Figueroa 
Street  Undercrossing  is  due  to  the  fact  that  arch  action  of  the 
deck  does  not  occur  because  the  amplitude  of  inward  motions  is  in- 
sufficient to  cause  closure  of  the  large  (6  in.)  joint  gap  in  this 
bridge  as  indicated  in  Figs.  74,  82  and  89. 

From  the  above  comparison,  it  appears  that  for  curved  bridges, 
a  bigger  joint  gap  is  desirable  in  order  to  avoid  biased  seismic 
responses  and  to  prevent  the  shake-down  situation  under  severe  shaking 
conditions . 

B.   FLEXURAL  YIELDING  OF  BRIDGE  COLUMNS 

For  both  5/14  South  Connector  Overcrossing  and  Figueroa  Street 
Undercrossing  Connector,  the  location  of  the  maximum  flexural  yielding 
in  the  columns  (Figs  103  and  104)  showed  that  under  moderate  shaking 
conditions  flexural  yielding  was  more  likely  to  occur  at  the  base  of 
the  columns  (Case  4  of  Fig.  103  and  Case  3  of  Fig.  104).   Under  severe 
shaking  conditions,  flexural  yield  hinges  also  developed  at  the  top 
of  the  prismatic  columns  (Cases  1,  2,  3,  and  5  of  Fig.  103  and  Cases  1 
and  2  of  Fig.  104) .   Columns  located  near  expansion  joints  were  shown 
to  be  more  likely  to  develop  yield  hinges  because  of  the  relatively 
higher  lateral  displacement  resulting  from  the  fundamental  transverse 
mode  of  vibration. 

The  damages  due  to  flexural  yielding  at  the  top  of  a  column  will 
penetrate  the  deck  and  hence  should  be  prevented.   An  effective  method 
of  achieving  this  is  by  flaring  the  upper  portion  of  the  column  as  in 
the  case  of  Figueroa  Street  Undercrossing  Connector. 

The  maximum  local  bending  ductility  factors  at  the  base  of  some 
columns  under  severe  shaking  substantially  exceeded  the  ductility  that 
could  actually  be  provided,  however,  under  moderate  shaking  they  re- 
mained within  tolerable  level  (  <5  )  which  could  be  accommodated 
by  well-designed  columns.   In  all  cases,  the  ductility  factors  required 
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by  shorter  columns  are  higher  than  for  the  longer  columns.   This 

results  from  the  fact  that  deck  segments  vibrate  laterally  almost 

as  rigid  bodies  causing  similar  lateral  displacements  of  all 

columns  and  hence  higher  ductility  requirements  of  the  shorter  columns. 

The  results  of  Case  1  (elastic  columns)  and  Case  2  (elasto- 
plastic  columns)  of  the  curved  Figueroa  Street  Undercrossing 
Connector  indicate  that  elastic  columns  substantially  increase  the 
acceleration  response  in  the  deck  without  significantly  affecting 
the  displacement  responses  of  the  system. 

C.   EXPANSION  JOINT  RESTRAINERS 

The  results  of  Cases  1,  2,  and  3  for  5/14  South  Connector 
Overcrossing  illustrate  the  effectiveness  of  longitudinal  restrainers 
in  reducing  the  seismic  response  of  the  bridge.   It  is  clear  from 
these  results  that  longitudinal  restrainers  at  expansion  joints  are 
essential  in  curved  bridges  with  several  expansion  joints.  The 
results  of  Case  2  illustrates  the  analytical  inadequacy  of  the  weak 
short  tie  bars,  a  fact  exemplified  by  the  bridge  during  the  San 
Fernando  earthquake.   Cases  1  and  5  with  strong  short  tie  bars  indicate 
sufficient  strengths  and  stiffnesses  but  inadequate  ductilities  of 
these  bars.  The  results  of  Table  27  showed  that  the  required  tie 
ductility  factors  of  the  strong  short  ties  greatly  exceeded  the  levels 
that  these  ties  can  accommodate  even  under  moderate  shaking  (Case  4 
in  Table  27) .   Hence  short  tie  bars  are  unsuitable  for  use  as  longitudinal 
restrainers  at  the  expansion  joints  of  long  curved  bridges  because  of 
insufficient  reserved  ductility. 

The  long  cable  longitudinal  restrainers  used  at  the  expansion 
joint  of  Figueroa  Street  Undercrossing  Connector  had  both  sufficient 
strength  and  ductility  to  resist  severe  shaking  conditions  (Cases  1, 
2  and  straight).  In  all  cases,  the  results  showed  that  these  cables 
did  not  experience  yielding  indicating  that  the  length  of  the  cables 
is  more  than  required  for  resisting  severe  shaking.  A  reduction  of 
2  in  length  would  have  still  enabled  the  cables  to  resist  severe 
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shaking  elastically. 

The  required  length  L  of  a  longitudinal  restrainer  tie  can 

be  estimated  from  the  expansion  joint  ledge  dimension  A  and  the 

Y 
yielding  strain  e  of  the  tie  by  the  relation 


LT=  — ?  (116> 

MT    T 


where  A  is  the  tie  gap  and  u  is  the  tie  ductility  factor  (Eq.  115)  which 
is  a  function  of  its  material  and  size.  A  can  be  estimated  from 
the  maximum  horizontal  displacement  response  of  the  bridge  and  the 
geometry  of  the  transverse  fundamental  mode  of  vibration  for  the  linear 
zero  friction  case. 

The  purpose  of  using  vertical  restrainers  at  expansion  joints  is 
to  prevent  spans  lifting  from  their  supports.  An  examination  of  the 
time  history  response  of  the  forces  in  the  vertical  restrainers  at 
expansion  joints  Nos.l  and  2  of  Cases  1  and  2  of  5/14  South  Connector 
Overcrossing  as  shown  in  Figs.  34  and  45,  indicates  that  under  severe 
shaking  tensile  forces  occur,  although  infrequently.   The  possibility 
of  tensile  forces  in  the  vertical  restrainers  indicates  the  desirability 
of  the  use  of  the  restrainers  to  prevent  vertical  separation  of  spans 
at  the  expansion  joint. 

D.   INFLUENCE  OF  DECK  CURVATURE 

The  influence  of  deck  curvature  on  the  dynamic  characteristics  of 
bridges  was  illustrated  by  the  mode  shapes  of  vibration  shown  in 
Figs.  68  and  93  for  the  curved  and  straight  Figueroa  Street  Under- 
crossing  Connector.   It  is  clear  from  the  results  presented  that 
horizontal  curvature  of  the  deck  causes  coupling  between  the  two 
horizontal  components ,  and  superelevation  of  the  deck  causes  coupling 
between  vertical  and  horizontal  components . 

The  effect  of  deck  curvature  on  the  seismic  response  of  the 
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bridges  was  shown  in  the  time  history  responses  presented  for  the 
straight  and  Case  2  of  the  curved  Figueroa  Street  Undercrossing 
Connector  (Figs  84-88  and  76-81) .   Ground  excitation  in  the  trans- 
verse direction  causes  greater  response  in  the  straight  bridge  which 
in  part  is  attributable  to  the  greater  lateral  stability  of  the 
curved  bridge  resulting  from  the  curved  distribution  of  supporting 
columns  and  abutments. 

E.  INFLUENCE  OF  VERTICAL  GROUND  ACCELERATIONS 

Vertical  ground  accelerations  have  an  important  effect  on  the 
vertical  response  of  all  bridges.  They  also  significantly  affect  the 
horizontal  response  of  curved  bridges  (Cases  1  and  5  of  5/14  South 
Connector  Overcrossing)  because  of  the  coupling  between  the  horizontal 
and  vertical  modes  of  vibration. 

Vertical  ground  accelerations  cause  large  vertical  oscillatory 
motions  in  the  deck  resulting  in  a  removal  of  the  vertical  compressive 
reaction  at  the  expansion  joints.  This  in  turn  removes  the  horizontal 
Coulomb  friction  force  which  affects  the  horizontal  resistance  of  the 
structure.  The  vertical  accelerations  also  change  the  axial  load  in 
the  columns  and  therefore  affect  the  ultimate  bending  strength  of 
the  columns  due  to  yielding  interaction  effect. 

F.  INFLUENCE  OF  COULOMB  FRICTION  AT  EXPANSION  JOINTS 

The  results  of  modal  analyses  presented  for  5/14  South  Connector 
Overcrossing  (Figs.  17  and  18)  clearly  showed  that  infinite  friction 
at  the  expansion  joints  substantially  increased  the  lateral  stiffnesses 
of  the  bridge  when  compared  to  the  case  with  zero  friction.   The  case 
of  zero  friction  is  considered  to  be  more  practical  and  realistic, 
because  the  Coulomb  friction  force  corresponding  to  the  friction  co- 
efficient normally  existing  at  expansion  joints  is  relatively  small  when 
compared  to  the  axial  force  which  would  result  from  the  arch  action  of 
the  deck  in  the  infinite  friction  case. 
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Because  Coulomb  friction  force  is  relatively  small  and  acts 
only  at  a  limited  nusttoer  of  expansion  joints,  its  influence  on 
the  total  bridge  response  is  small. 

G.   INFLUENCE  OF  OVERALL  STRUCTURAL  ARRANGEMENT 

Although  longitudinal  restrainer  ties  are  effective  in  reducing 
the  amplitude  of  seismic  response,  they  do  not  affect  the  dynamic 
characteristics  of  the  bridge.   This  was  illustrated  by  the 
essential  similarity  between  the  shapes  of  time  history  response  of 
Case  1  and  2  of  5/14  South  Connector  Overcrossing.  Longitudinal 
restrainers  are  secondary  structural  elements  which  are  effective 
only  when  the  primary  horizontal  seismic  response  is  sufficiently 
large  to  activate  them. 

From  the  results  presented  for  5/14  South  Connector  Overcrossing 
and  the  curved  Figueroa  Undercrossing  Connector,  it  is  concluded 
that  the  major  influence  on  the  seismic  response  characteristics  is  the 
arrangement  of  supporting  columns  of  each  deck  segment.   For  sufficient 
earthquake  resistance  of  a  bridge  structure,  each  deck  segment  should 
be  supported  by  at  least  2  columns  so  that  its  lateral  stiffness  is 
reasonably  adequate,  and  that  stable  oscillatory  response  of  the  structure 
can  be  achieved. 
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IX   CONCLUSIONS  AND  RECOMMENDATIONS 


A.   CONCLUSIONS 

Based  on  the  analytical  results  presented  for  the  bridge  structures 
sttidied  in  this  investigation,  the  following  conclusions  may  be  deduced: 

(1)  The  basic  mathematical  models  and  the  analytical  procedures 
as  well  as  computer  programs  developed  in  this  investigation 
provide  a  rational  and  effective  means  for  determining  the 
seismic  response  of  curved  and  straight,  multiple-span,  long, 
highway  overcros sings  of  the  type  considered  herein. 

(2)  The  collapse  of  5/14  South  Connector  Overcrossing  during 
the  San  Fernando  earthquake  was  due  to  the  inability  of 
the  longitudinal  restrainer  ties  to  constrain  the  joint 
separation  at  expansion  joints  Nos .  1  and  2  thus  causing 
the  central  spans  between  these  expansion  joints  to  fall  off 
their  supports  resulting  in  the  collapse  of  the  spans  along 
with  column  No.  4. 

(3)  The  proposed  curved  Figueroa  Street  Undercrossing  Connector 
can  safely  withstand  severe  ground  shaking  of  an  intensity 
similar  to  those  experienced  during  the  San  Fernando  earthquake* 
however,  some  flexural  damage  is  likely  to  occur  at  the  bases 

of  the  columns  supporting  the  central  spans. 

(4)  The  strength  and  distribution  of  bridge  columns  in  each  deck 
segment  between  expansion  joints  are  the  primary  factors  which 
most  significantly  affect  the  seismic  response  characteristics 
of  the  bridge  structure.   Thus  expansion  joints  in  the  deck 
should  be  located  in  a  manner  such  that  each  substructural 
segment  posses  a  reasonably  adequate  amount  of  lateral  stiff- 
nesses . 

(5)  Under  severe  shaking  conditions  similar  to  those  experienced 
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during  the  San  Fernando  earthquake ,  the  maximum  horizontal 
acceleration  response  at  the  deck  level  substantially  ex- 
ceeded those  values  specified  by  the  "Bridge  Planning  and 
Design  Manual",  California  State  Division  of  Highways  [5j 
which  was  in  effect  prior  to  the  San  Fernando  earthquake 
of  February  9,  1971. 

(6)  Because  of  their  effectiveness  in  reducing  the  amplitude 
of  the  seismic  displacement  response  of  curved  bridges, 
longitudinal  restrainer  ties  should  be  provided  at  expansion 
joints.  They  should  have  sufficiently  high  strength  and 
ductility  to  constrain  the  joint  separation  within  the  limit 
of  the  supporting  ledge  dimension.   Vertical  restrainers 
should  also  be  provided  at  expansion  joints  in  order  to  pre- 
vent the  lifting  of  spans  from  their  supports. 

(7)  For  curved  bridges,  a  bigger  expansion  joint  gap  appears  to 
be  desirable  to  prevent  biased  seismic  response  and  to  avoid 
a  shake-down  situation  under  severe  shaking  conditions . 

(8)  The  required  expansion  joint  ledge  dimension  and  joint  gap 
can  be  estimated  from  the  maximum  horizontal  displacement 
obtained  from  a  seismic  analysis  and  the  geometry  of  the 
transverse  fundamental  mode  of  vibration. 

(9)  Under  moderate  to  severe  shaking,  the  overall  seismic  response 
of  a  bridge  structure  is  considerably  affected  by  vertical 
ground  accelerations.   The  effect  is  particularly  noticeable 
for  curved  bridges  in  which  coupling  exists  between  vertical 
and  horizontal  modes  of  vibration.   Thus  the  vertical  component 
of  ground  acceleration  should  always  be  included  in  a  seismic 
analysis. 

(10)  For  a  normally  proportioned  bridge  structure  under  moderate 
to  severe  shaking ,  flexural  yielding  is  more  likely  to  occur 
at  the  base  of  the  columns.   Yielding  in  columns  should  be 
avoided  under  moderate  shaking  but  can  be  permitted  to 
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limited  levels  under  severe  shaking  conditions.   It  is 
therefore  important  that  column  reinforcement  should  be 
sized  and  detailed  in  such  a  manner  that  the  columns  can 
develop  full  yield  strengths  and  maintain  their  yield 
capacities  throughout  their  specified  ranges  of  ductility. 
Flexural  yielding  should  be  prevented  from  developing  at 
the  top  of  columns  in  order  to  avoid  damage  penetrating 
the  deck. 

(11)  The  effect  of  Coulomb  friction  forces  on  the  seismic  response  of 
the  type  of  bridges  considered  is  small,  for  the  friction 
coefficient  values  normally  existing  at  the  expansion  joints . 
Thus  for  linear  analysis,  a  mathematical  model  assuming  zero 
friction  at  expansion  joints  can  be  used  as  a  realistic 

linear  bridge  model. 

(12)  Linear  seismic  response  analysis  provides  a  reasonable 
estimate  of  the  maximum  displacement  response  of  a  bridge 
system;  however,  substantial  error  may  result  in  predicting 
the  internal  forces  in  the  structure.  Thus  it  is  essential 
that  nonlinear  seismic  analysis  be  carried  out  for  major 
bridge  structures  under  severe  shaking  conditions. 

B.   RECOMMENDATIONS 

Due  to  the  unique  dynamic  characteristics  of  bridge  structures 
considered  in  this  investigation,  it  is  recommended  that  major  highway 
overcrossings  should  be  designed  on  the  basis  of  dynamic  analyses 
using  appropriate  nonlinear  mathematical  models  to  represent  the 
structural  systems.  The  dynamic  analysis  should  include  both 
vertical  and  horizontal  components  of  ground  motions  suitable  to 
the  respective  site.  In  determining  the  seismic  response  of  long  bridges, 
the  "out-of-phasing"  of  the  multiple  support  input  ground  motions  may,  in 
some  cases,  be  significant.  Therefore,  every  effort  should  be  made  to 
measure  spatial  correlations  as  well  as  time  correlations  for  strong 
ground  motions. 
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